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ABSTRACT 



In this research a general purpose digitally controlled 
analog filter is presented. The novel design is a cascade of 
second order sections that are individually programmed to 
achieve any filtering topologies. Two-binary words are used 
to control the pole frequency co p and selectivity Qp of each 
section independently. Each second-order section 
is a Generalized-Immittance Converter (GIC) biquads which are 
known for their high stability and low active and passive 
sensitivity. CMOS switches are used to electronically 
relocate the minimvim number of passive elements to achieve 
function programmability. Switches are also used to select 
the number of cascaded sections to realize higher order 
transfer functions. 
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I. 



INTRODUCTION 



A. THE NEED FOR AN ACTIVE PROGRAMMABLE FILTER 

The availability of an analog filter with digitally 
controlled "programmable" coefficients has been the goal of 
many researchers due to its several attractions. One 
possibility of a compact, versatile analog filter under 
remote control opens up many novel and independent 
application areas. Also, when a programmable filter is 
combined with a permanent referenced memory which is 
user-programmable, this would form an economical and 
versatile device for dedicated stand-alone applications. The 
need for such a device was motivated by advancement in thick 
and thin film technologies and continuous upgrading of 
systems specifications to take advantage of the available 
technologies to the limits. 

Linear analog filtering finds many applications, such as 
speech processing (recognition or synthesis), geology, 
instrumentation, communications, process control, adaptive 
balancing, etc. There has been much emphasis on performing 
the filtering function digitally, largely because of the ease 
of varying and optimizing the transfer function. However, 
for many reasons, such as cost, size, signal processing 
complexity, and bandwidth, it would be desirable to perform 
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the filter function with linear components yet retain the 
flexibility of varying the filter parameters digitally. 

Recently, the advantages of combining linear components 
(operational amplifiers, resistors and capacitors) and non 
linear elements (switches) have been demonstrated using 
switched capacitor techniques [36-38]. In this research, we 
are presenting the results of realizing a continuous active 
device using linear elements and switches controlled by 
digital signals to achieve a fully programmable filter [32]. 
Several programming features of the proposed filter are 
reported. The first feature is the ability of the networJc to 
realize the most common filtering functions (function 
programmability) namely: Low Pass (LP) , Band Pass (BP) , High 
Pass (HP), All Pass (AP) and Notch (N) functions, using the 
minimal set of elements. The second feature is the ability 
of the networJc to program (independently) the )cey parameters 
of the filtering function chosen (parameter programmability) 
namely: the pole resonant frequency (u^p) and selectivity 
(Qp). Finally, the ability to program the networ)< to cascade 
several sections to achieve higher order filer. All of the 
above programmability features are performed independently to 
realize a universal filtering networl<. In order to 
demonstrate the idea of this research, it is necessary to 
introduce some theoretical bade ground. 
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B . LINEAR SYSTEMS 



The box at Fig. 


(1.1) 


illustrates 


the 


concept of 


the 


linear system. In 


the 


time domain, 


the system 


is 


characterized by its 


impulse response 


y(t) 


, which is 


the 


output signal y(t) produced 


in response 


to 


an impulse . 


for 



an arbitrary input signal u(t), the output signal y(t) is 
given by the well known convolution integral. 






aLL-l) XLO 
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vj 
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ULL) 
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» !•! “ Black Box Concept of Linear System 

The system transfer function G(s) is the Laplace 
transform of g(t), thus: 



CsL'j'\ ~ tj LL'i (ixp C-b k.') 
iwlntif e. S = iT \ cjo 



( 1 . 2 ) 
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Under the laplace transformation 1.1) becomes; 



= dts:) • \Jc5D 



(1.3) 



so that the transfer function G(s) is the ratio of the output 
variable to the input variable. 



' The most general types of linear systems are consisting 
of a finite number of lumped, linear and time-invariant 
elements. The system is characterized by an Nth order 
ordinary linear differential equation which results, in most 
general cases, a transfer function G(s) which is a real 
rational polynomial function of the complex variable S. Thus 
we can write G(s) in general as 



where Pk and Ik are real numbers so that G(s) is real for 
real s, and the roots of the polynomial P(s) and E(s) must be 




(1.4) 
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real or occur in conjugate pair. 

Also by proper multiplications G(s) can take the form 

TV CS-52k) 

( 1 . 6 ) 

TT (S-SnO 

K 

where C is a constant extracted from G(s) such that E(s) and 
P(s) become monic polynomials (leading coefficients equal to 
unity) , Szk are the transmission zeros and Snk are the 
natural modes at the system. 

C. FILTERS AS A SPECIAL CLASS OF LINEAR SYSTEMS 

Linear systems can be distinguished into "SPECTRAL 
SHAPING NETWORKS" and "FILTERS." The role of filters is one 
of selecting signals while the role of spectral shaping 
networks is that of modifying the input signal spectrum in an 
arbitrary, but predescribed manner. Specifically, we desire 
that a filter should do as little as possible shaping on 
signals in its passband; any shaping it is considered a 
distortion of the signal. On the other hand, networks which 
perform pulse forming fall within the spectral shaping 
category. 

D. ACTIVE FILTER FUNDAMENTALS 

The distinction between passive and active filters is 
that the first do not require a power source to perform their 
function while the second do. The motivation behind active 
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RC filters lies in the desire to have inductorless filter 



realizations. It is well known that of the three passive 
R,C&L elements, the inductor is the most nonideal one. This 
is especially true at low frequencies, where inductors become 
quite bully and have increased losses or equivalently lower 
Q-f actors . 

E. GENERALIZED IMMITTANCE CONVERTOR 

One of the methods of active RC filters design consists 
of simulating the inductances in the LC ladder by active RC 
networks. This simulation can be based on the principle that 
we want to find a one port network having an input impedance. 

Zll = s*L 

Various active elements as well as synthesis procedures 
employing them have been proposed [1-9]. A partial list 
includes : 



(1) 


Negative Impedance Converter 


(NIC) 


(2) 


Negative Impedance Inverter 


(NIV) 


(3) 


Postive Impedance Converter 


(PIC) 


(4) 


Gyrator 




(5) 


Generalized Impedance Converter (GIC) 


(6) 


Curent Conveyor ( CC ) , and 




(7) 


Operational Amplifier (OA) 





Although the introduction of these elements has 
stimulated research in the area of active network theory, 
very few elements have made their way to large scale 
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production to become available as off-the-shelf items. The 
reason for this is mostly an economic one. For a device to 
become available at low cost, it has to be used in 
substantially large quantities. It follows that such a 
device has to be versatile enough to be of use in a number of 
applications, of which active filter design is only one. The 
analog circuit design area has found these attributes in the 
IC operational amplifiers (Op. Amp.) 

The IC Op. Amp. is currently the most popular linear 
active element. It is available from a large number of 
manufacturers, at reasonable cost and with good performance 
characteristics. Furthermore, elements engineers have become 
accustomed to the use of Op. Amp. It is therefore, only 
natural that the Op. Amp. is becoming the most popular active 
element in the design of active RC filters, and can be found 
in NIC's, PIC's, QIC's, and other circuit realizations. 
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II. THEORETICAL ANALYSIS 



A. BIQUADRATIC TRANSFER FUNCTIONS 

The filter as a special class of linear system has a 
transfer function expressed in a polynomial quotient form 
given as 



where Pk and gk are real numbers to that T(s) is real for s, 
and the roots of the polynomials P(s) and E(s) must either be 
real or occur in conjugate pairs. Also, in general, the 
degree of the numerator (deg.[P(s)], = M) is less than or 
equal to the degree of denominator (deg[E(s)] = N) and the 
roots of E(s) are in the open-half S-plane. The E(s) is 
known as the characteristic polynomial or natural mode 
polynomial of the linear system, and the degree of E(s), that 
is N, is the order or degree of the system. 

A general second-order transfer function or "biquad" 
function may be written as 




( 2 . 1 ) 
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^ ^ S 1- Fcs^ 



tc>:> 



( 2 . 2 ) 



where P(s) is the loss-pole, or more appropriately here, the 
transmission-zero polynomial, and E(s) is the natural pole 
polynomial mode as discussed above. It is a usual practice 
to express the denominator in terms of _p and Qp, where _p is 
the natural-mode or resonance frequency and Qp is the 
natural-mode or quality factor. Thus (2.2) becomes 



The numerator coefficients determine the location of the 
transmission zeros and hence, the type of filter function the 
biquad provides. Special cases of interest are; 

1 . Low Pass (LP) 



) _ T; S f/^S 1- to 




(2.3) 



For which Pl=P2=0, thus two transmission zeros are at 



infinity 




(2.4) 
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2 . High Pass (HP) 



For which p0=pl=0, thus two transmission zeros are at 
infinity, and 



\Cs^~ 






UOP y 

1 S V 00^ 

Qp 



(2.5) 



3 . Band-Pass (BP) 

For which Po=P2=0,thus one transmission zero is at 
infinity while the other is at the origin, and 









5^ 



s -t- U) P 

cpp ' 



( 2 . 6 ) 



4 . Notch (N) 

For which P1=0 and the two transmission zeros are at 
5 rl jo;,! ^ 10 ^/^ (Up (depending if we have low-pass-notch 

or high-pass-notch) , leading to 



I CS^ : I 









(2.7) 



Qp 
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5. All Pass (AP) 

For which the pair of zeros are at the mirror image 
location of the pair of poles, that is 



B. SENSITIVITY FUNCTIONS 

A concern about the design of a filter is how close the 
resulting response will be to the ideal or desired function. 
The reason for response deviation from the ideal is the 
finite tolerances of the RC elements, as well as the nonideal 
performance of the active elements. In the latter case, not 
only the gain changes or tolerances have to be considered but 
the effect of the "limited amplifier bandwidth" on the filter 
response must also be evaluated. Although effects of initial 
component tolerances may be "trimmed out" during the initial 
filter alignments or tuning process, a sensitive design will 
deviate from the required specifications as time process, due 
to component variations with temperature, aging, humidity, 
etc. Note also that a sensitive design might be extremely 
difficult to tune in the first place, or the initial 
adjustment will be quite uneconomical. 

The answer to tolerance question can be obtained through 
sensitivity studies. Considerable emphasis has been placed, 
in the active filter literature, on the study of sensitivity 




( 2 . 8 ) 
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functions and relations. some of the most useful and widely 
accepted sensitivity functions are the: 

1 . Magnitude Function Deviation 

Assure a filter designed to meet a certain magnitude 
characteristics T(s) or T(j>^). One is concerned with the 
deviation in jT(jO)j, that is/l| T( jv^) | both in passband and in 
stopband. Usually it is desirable to express the expected 
deviation in dB. The deviation D(uJ>) dB in the magnitude 
function may be evaluated as follows. Let the function iT(j^».i 
change to [|T(j-xi)| -t jT(j )|], then 



Ddtu) r: QOlorj 1 ^ ^ ^ 


(2.9) 


1 1 Cj uD^ 1 





or 



Dlid) - b bZ ^VL 



4 . + 



L 



i ^U^)l 



,( <4b) 



( 2 . 10 ) 



and for small variability (2.10) can be approximated as, 



Pluo) 2^ 8.68 ^ I ^ U 

I 1 






( 2 . 11 ) 
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Thus, the deviation in the magnitude response in nepers 
is equal to the per unit variability in the magnitude of the 
transfer function. The problem now reduces to that of 
evaluating the per unit change in jT(jvP)|. This is not an 
easy problem since T(juJ) is a function of many elements with 
different tolerances and tolerance statistics. Furthermore, 
the per unit change is function of frequency. 

2 . Classical Sensitivity 

Lets recall the definition of the classical 
sensitivity, where y is a variable of interest, usually a 
function of many parameters of which x is one, then 



5 ," ^ ^ :L 



( 2 . 12 ) 



Note that from the above definition, is the limit to as 
Dx -->0. thus, for small variations. 




(2.13) 



The usefulness of the classical sensitivity function 
is evident from (2.13). The per unit or percentage change in 
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y, due to a given per unit or percentage change in x, can be 



easily obtained by multiplication with S, i.e.. 




(2.14) 



3 . Gain Sensitivity Product 

An important consideration in the evaluation of the 
sensitivity of a filter parameter as considered in [2] with 
respect to the closed loop gain is the tolerance on the 
closed loop due to the open loop gain variability. 

Thus, the gain-sensitivity-product, G.S^ is defined 

as: 



6. S K 5 

1^ K 



(2.15) 



where k is the closed loop gain. We can extend (2.15) to the 
open loop gain as: 



(s S'^ Ao 

Ao '^0 



(2.16) 



and also we can note that: 
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(2.17) 



.g.s;=gs 



Ao 



and thus then ultimate good is the variability or tolerance 
rather than the sensitivity, the gain-sensitivity product is 
a better index for comparing different designs. 

4 . Determining the Variability of the Transfer Function 
Amplitude 

Assure that the active filter has t resistors, 
tTL capacitors and n amplifiers. Let the amplifier k have an 
open loop gain Aok and, possibly, a closed loop gain Kk. The 
variability of the magnitude function as given by the 
Reference [10] is: 



1 T(_jud)1 




m. 







s 










+- 



A: 

K--1 






(2.18) 



Note that each of the sensitivity functions is (2.13) 
is a function of frequency. In a high order filter 
realization, the different sensitivity functions might be 
difficult to evaluate. 
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C. PROPOSED QIC FILTER ANALYSIS 

In order to obtain the transfer functions of the proposed 
programmable filter [13] shown at Fig. (2.1), nodal analysis 
was used as follows: 

(1) The circuit of Fig. (2.1) was replaced by the one of 
Fig. (2.2) in which the two operational 

amplifiers were replaced by the two equivalent dependent 
voltage sources as can be seen, every element (node, 
admittance, voltage source, etc.) was labeled and every 
element was associated to a current direction and a voltage 
polarity. 

(2) The kirchoff current low was written for every node 
except: 

(a) The reference, 

(b) any node connected to the reference by a voltage 
source . 

Node K.C.L. 

1 i7=-i3 (2.19) 

2 i2=i5+i6 (2.20) 

3 

4 

5 i4=i7+i8 (2.21) 

6 
7 
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Fis. 2.1 



Active Filter Con f igurat ion Using the GIC 
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:3 




Fig. 2.2 - The Circuit of 2.1 prepared for Nodal Analysis 
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( 3 ) Every admittance current was expressed in terms of 
nodal voltages; 



il = Y1[V3 “t/i] 


(2.22) 


i2 = Y2h/3-V2] 


(2.23) 


i3 = Y3[o/4-Vl] 


(2.24) 


i4 = Y4[V4-VS] 


(2.25) 


i5 = Y5 [t/2-1^] 


(2.26) 


i6 = Y6[V2] 


(2.27) 


i7 =[Y7[^^5-V6] 


(2.28) 


i8 =Y8[lJ5] 


(2.29) 



( 4 ) The source dependencies were listed expressed in terms 
of nodal voltages, 






and then (substitute) where necessary. 

(5) A matrix equation having the unknown voltages V3, V4, 
and V5 related with the desired transfer functions Tl, T2, 
T3, of the filter was obtained as: 



■^1- f' 1 t — 

L J ^2 










a33 

1 






'^1 1 


Ai A, Ai 






i4 










'{ f t)., 



(2.32) 



0 
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and by substituting 






T3 = '“j- 



^\)4 

and ' 3 - 

"OiiO, "U 

the above matrix equation takes the following form; 




'<5 





,H_5 '(£_ 

h\ f\\ l\\ 



0 



'^4 









t) 






.-1 






! -^5 



-4 






(2.33) 



So the above matrix equation gives the 3 different 
responses of the programmable filter as functions of the Yl, 
. . . Y7 admittances and Al, A2 (Gains) the two Op. Amps. In 
ideal case [Al=A2=oo] the equation takes the form: 



Hi 



X5 



-IvHi 







*1 

•s 



0 






T:? 




0 



X4 






'5 



—I 




(2.34) 



which expressed the different realizations discussed by [31]. 
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Now we can express the above matrix equation for the 
nonideal case in which A1 and A2 are finite and frequency- 
dependent I where W1 and W2 are the Gain 

Bandwidth Products (GBWP) of the Op. Amps. Al, A2, 

respectively. 





. ^-4 s 




-is 


i MCs'j . 

1 


1 ^ 




,1 ['Isv'Cc] 

vMj 


Wi - 


1 

>i 'i. \l 

- -I- '■•3 '4 


T , - 'J' 


0 






-\<4 
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1 . Low Pass (LP) Realization 

Can be obtained through T2(s) by substituting 
Y6=Y7=0, and if we consider the ideal case where (Al=A2=oo) , 
then 
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(2.36) 



Hi ^5 1- '{ft 



If we further substitute the values of the remaining 
admittances as proposed in [4], that is Yl=Gl, Y2=sC2, 
Y3=sC3+G3, Y4=G4, Y5=G5 , Y8=G8 then (2.36) takes the form 



(2.37) 

(Sv(r4 S>^ 



which is the form of a low pass transfer function as 
indicated by (2.4) where 
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(2.38) 



_ <^, &A 6s 
6z}C:^ 



(2.39) 



and 






(2.40) 



Op C3 



2 . High Pass (HP) Realization 

This one is obtained through T1 ( s ) substituting the 
following values of the admittances: Y1=G1, Y2=G2, Y3=sC3 , 
Yr=G4 , Y5=0 , Y6=G6, Y7=sC7, and Y8=G8, and if we assume ideal 
case (Al=A2-> 00), the 






(2.41) 






expresses the H.P. filter of (2.5) where 




(2.42) 



Qb 



(2.43) 



<Op c> 
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and 




(2.44) 



3 . Band Pass (BP) Realization 

This is also derived from Tl(s) by substituting the 
following values of admittances: Y1=G1, Y2=G2, Y3=sC3, 

Y4=G4, Y5=0, Y6=G6, Y7=G7, and Y8=sC8, and if we assume ideal 
case (A1=A2 — >oo), then 



which is the form of a BP transfer function as given by (2.5) 
where 




(2.45) 







(2.46) 









(2.47) 







(2.48) 



Cq (?2 
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4 . Notch (N) Realization 



T2(s) expressed the Notch response if the following 
substitutions have been made. Y1=G1, Y2=G2, Y3=sC3, Y4=G4, 
Y6=0, Y5=G5, Y7=sC7, Y8=G8 , and again assuming ideal case 

(A1=A2 — >oo) , 



^ Q-2 +■ S <5-4-G5 



(2.49) 



which is the form of a N transfer function as defined by 
(2.7) where 



11)^ - 



6-5 (2<S-4 
^3 Cl Gz. 



(2.50) 



(C\ <$~4 Gs 
Ci>CiQ2, 



Qp C7- 



(2.51) 



(2.52) 



5 . All Pass (AP) Realization 

This is derived from Tl(s) with the following 
admittances substitutions. Y1=G1, Y2=G2, Y3=sC3, Y4=G4, 
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Y5=G5, Y6=0, Y7=sC7 and Y8=G8 , and if once more we assiime 
ideal case, then 



^ sCs Gs G5 

^3 C'}-G'2+ S (^3(32 G© 1 - G"!? (Ti G 4 



(2.53) 



which is the response of an All Pass filter as (2.8) 
indicates where: 



(iO? 

Qp 



($e 

Cj 



C ^O-C 



noo>T)i niiDoo ^hose ^ 



(2.54) 



0\->6 

C3 Ct- 



) 



(2.55) 



Table 2.1 shows all the realizations proposed by 
[31]. In our research for designing a programmable filter 
the No. 1, 3, 7, 9, and 12 realizations were used since they 
offer the minimum admittance elements . change to shift from 
one to another. 

D. SENSITIVITY ANALYSIS 

Consider the CGIC circuit shown in Figure (2.2(a)) [12]. 
Assuming ideal Op. Amps., the chain matrix of the CGIC can be 
obtained as 
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(o) 



CGIC I \ 0“^ 








Fig. 2.3CoOThe CGIC Implimentation Using 
Op, Amps. 

(b) Symbolic Representation of the 
CGIC with created ports 3G and 4G 

(c) The Basic Conf irguration . 
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2.1 - Elements identification for 

the Most Commonly Used Transfer Fuctions 
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1 



0 



L^l = 



(2.56) 



0 WCS) 



where h(s), the admittance conversion function, is given by: 



Two new ports can be created across terminals 3G and 4G 
as shown symbolically in Fig. (2.3 (b)). 

A synthesis procedure is now described which uses the 
configuration of Fig. (2.3(c)). The transfer functions 
between the input and output terminals 2, 3, and 4 are 

readily obtained as 

"OihK -- r (2.58) 



il ii l'(\ U 



(2.57) 



r 1 '(si. U'J'oh.A') V- Oci) 



(2.59) 



T%r H^'slDo:) 



(2.60) 
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where 

~ Ko V Ng V kcs)C 



(2.61) 



The conversion function h(s) and Y5-Y8 can be selected in 
many different ways and it is found that any second-order 
transfer function can be realized [12]. 

Letting 

~ sC 2 . 1 - Qt. 

where i=l, 2 , 3, or 4, we have from (2.57) 

lie's) r C sCl ^63) / CsCi '-Gx ) ()sC4 + G4)) ( 2 . 62 ) 



Clearly by omitting one or more conductances and/or 
capacitances a number of specific conversion functions can be 
generated. 

Most frequently, filters are designed by using 
Butterworth, Chebychev, Bessel or elliptic approximations in 
which the transmission zeros are located at the origin, 
imaginary axis or at infinity. Consequently, the. transfer 
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function can be expressed as a product of a second-order 
transfer functions of the form 



1 cs) r 



O2 Oh s + ao 
bj S+fc>H^+-bo 



(2.63) 



where al=a2=0 , ao=al=0, ao=a2=0 or al=0 for low pass (LP), 
high pass (HP), band pass (BP) or notch (N) section, 
respectively. 

The coefficients ai's of T(s) for these sections are all 

positive. These sections can be realized by choosing h(s) in 

2 

a simple manner such as k^s, k 2 S, k^s +k^s or their 
reciprocals. The different ki's (i=l, 2, 3, or 4) are 
positive constants. By comparing (2.58-2.61) with (2.63), 
circuits 1-10 in Table 2.1 can be obtained. Circuits 3, 4 
and 7 can be regarded as realizations of simple RLC networks 
[31]. 

All pass transfer functions are often needed for delay 
equalization and these can be realized by using second-order 
transfer functions of the type give by (2.63) where a2=b2, 
al=bl, and ao-bo. Second-order sections of this class can be 
obtained from circuits 11 and 12 of Table 2; 1. 

Figures (2.3(a)), (2.3(b)), and Table 2.1 show that with 
the exception of circuit 10, the response is obtained from 
the output of an operational amplifier. Owing to the low 
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output resistance of the amplifier, any number of sections 
can be cascaded without isolating amplifiers. 

An important criterion of a realization is its 
sensitivity to element variations. The pole Q factor and the 
undamped frequency of oscillation from the transfer function 
of (2.63) are defined as 



For a Notch section, the Notch frequency is defined by 





(2.64) 




(2.65) 



and the multiplier constant can be talcen to be 




( 2 . 66 ) 



for (tn>ii>p or tt>n<iiip, respectively. 



Similarly for the LP, HP, and BP sections 




(2.67) 
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For an All Pass Section, let 






Q'i - ^Qo ^ (J03=(^Q o I Q 2^ 

O \ 



'U 



Ood 



( 2 . 68 ) 



The sensitivity of a quantity x with respect to 
variations in an element e is given by 

For ideal amplifiers, the use of ( 2 . 64 ) - ( 2 . 68 ) and Table 
(2.1) leads to 







( 2 . 69 ) 



where x represents any one of the quantities defined by 
( 2 . 64 ) - ( 2 . 68 ) and e represents any capacitance of 
conductance. In addition, it can be shown that 
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and 




(2.70) 



For amplifiers with a finite open-loop gain A, according 
to [31] the circuit of (Fig. 1(c)) gives 





(2.71) 



where k=2 , 3 , 4 , and 

=fi + r-i)(uv'2')(.hlAM- fi/Ai +Xl/A^^2 ^''3 '3,A2) I (2.72) 

Ti - " 

Consider realizations in which h(s) = kls, such as the 

circuits 7 (BP), 3 (HP), 9, 10 (N), where 

1'.-, (s'S ) S (s3>, X-i - ^4 ) ^-5 , \c Gb ) 1?- SC? Ovii 

td-sCb -Gb 



(2.73) 



For real amplifier gains such that 

A\ ^kir^ko Ko y) L 



L^i^Xgives ^ 

Ub^) r Ft Hi V vHaVAo 
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(2.74) 



From (2.64), (2.73), and (2.74), the Q-f actor and the 



undamped frequency of oscillation can be obtained as 



Qp , qA \ ^ i \ / [ i 






C:3 A3 Ad 



j 




^ X5 ^ 

(<^T4(fs)Ao ' 



(2.75) 



where 

Op r \ (jJ5 *-0) CC I f-Cfc) (si64 j O? ^ 

dip- f ((y5^0) I CCit-C‘b)($:?C3 ^ 

'Ai.- 

/ 

^ -2 ' ((rl } / (5<t , 



(2.76) 



'^■'5- (fCi -(-(re^/6'f 

y^S ~ ood ^4-C^rAy^^-<4- 

The sensitivities of Qpa and pa with respect to the 
amplification Ao can be written as 




-t ]' X4 
.3^o L ^ '^'3 






li)V4 

Sao 



\ 

3Ao 





(2.77) 



The use of (2.77) and (2.78) leads to 




(2.78) 



by assuming that 



0 ^ ^'2Q^ 'KjjO ^ 

’ ')k( 



(2.79 



Eq.(2.78) reduces to 



"^Ao ' 




(2.80) 



Straight forward differentiation shows that 
when 



sQpa 

1 .) 



is minimum 







(2.81) 



From (2.79) and (2.81) the analysis is valid provided 
that >/ 1, which is clearly satisfied in 

practice. 

From (2.80) and (2.81) the minimum sensitivity to 
variations in Ao is derived as 






/Ao 



(2.82) 
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(2.83) 



The corresponding value of 



is given by 

Ao 



^Ao 



J 

Ao^p 



LP realizations as circuits 1 and 2 can be obtained by using 
a conversion function of the form or its reciprocals. The 
admittances Y1 to Y8 are chosen as 



in order to obtain a conversion function of the form 



Now Qpa and c^pa are obtained as 



..'h 



^ ‘ ' 

^ r /^aC? AoJ 



(2.85) 



i!.7 






( 2 . 86 ) 



Ih 



0..: f j (arr3'(^) ]■'" 



li' 



I (C-I 6-4 65' ; j (QC3 % )]' 






(2.87) 
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The sensitivity of Qpa with respect to variations in Ao 
can be minimized following the approach used earlier. It is 
found that for minimxjm sensitivitv C 31 ] 

^ ( 2 . 88 ) 

The minimum value of can be shown to be 

Ao 






/Ao 



(2.89) 



and the corresponding value of is given by 

Ao 






(2.90) 



The above sensitivity analysis can be extended to 
realizations using any other type of conversion function 
[31]. 

Equation (2.69) shows that the sensitivities to passive 
element variations are independent of the selectivity. 
Furthermore, the sensitivities with respect to variations in 
the amplifier gain are low. The proposed realizations are 
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seen to have similar sensitivity properties as the low 
sensitivity realizations reported in [21-27]. 

E. STABILITY 

It has been shown elsewhere [20] that some networks using 
QIC's can be conditionally stable where a circuit can lock in 
an unstable mode during activation (just after switching on 
the power supply) . In this section the stability properties 
of the configuration show in Fig. (2.2(c)) are examined. 

The natural frequencies of the circuit in Fig. (2.3(c)) 
are the zeros of the characteristic polynomial D9s) as given 
by (2.72). The differential open-loop gain of a frequency 
compensated Op. Amp., in a bounded frequency range i 

Ao JJl) c 1 C S ^ Jj c, ^ 



where Ao and wc are the d.c. gain and the cutoff frequency 
respectively, and 0 ^ Ao ^ /V.v'.o;, . In the frequency range 

'1)<<'-22c the amplifier gains A1 and A2 can be assumed to be 
real. For any second-order transfer function the coefficients 
of D(s) are seen to remain positive for any attainable pair 
of A1 and A2. this is due to the absence of negative terms 
in D(s). Therefore, the zeros of D(s) will remain the 

left-half s-plane and low frequency unstable modes cannot 
arise during activation. 
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III. PROGRAMMABLE GIC FILTER 



A. GENERAL 

Signal processing devices evolved considerably over the 
last several years. The progress was motivated by the 
advancement in film and semiconductor technologies, as well 
as the continuous upgrading of systems specifications to take 
advantage of the available technologies to the limits. 

Linear filtering finds many applications, such as speech 
processing (recognition or synthesis), geology, 
instrumentation, communications, process control, adapting 
balancing, etc. There has been much emphasis on performing 
the filter function digitally, largely because of the ease of 
varying and optimizing the transfer. However, and for many 
reasons, such as cast size signal processing complexity, and 
bandwidth, it would be desirable to perform the filter 
function with linear components, yet retain the flexibility 
of varying the filter parameters digitally. 

Recently, several advantages of combining linear 
components (amplifiers and capacitors) and nonlinear elements 
(switches) have been demonstrated using MOS switched 
capacitor techniques [31, 32]. Here, we are presenting the 
results of realizing a continuous active device using linear 
elements and switches controlled by digital signals to 
achieve fully programmable filters. 
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Our research addressed two different aspects of 
progranunability namely; 

(1) Programming the filter topology using a minimal set of 
elements to obtain any type of filtering function desired, 
e . g . , LP , HP , BP , N and AP . 

(2) Programming the filter's transfer function parameters, 
(pole resonant frequency iup and quality factor Qp) for a 
chosen type of filtering function. 

B. THE PROPOSED GIG PROGRAMMABLE FILTER 

The basic active network considered as the heart of the 
GIG programmable filter is the GIG structure [31] of 
Fig. (3.1), whose superior performance was established in the 
literature [10,32]. The filter transfer function was derived 
using loop analysis in Ghapter II. 

Table (3.1), illustrates that for any of the LP, HP, BP, 
N and AP realizations, five resistors, two capacitors and two 
Op. Amps., are required. also, the transfer function of each 
realization is shown. The passive elements are connected to 
the different nodes, shown in fig. (3.2), for the different 
realizations. A set of MOS bilateral switches controlled by 
a digital binary word, are used to interchange the elements 
to achieve the different types of filter realization shown in 
Fig. (3.3). The truth table of the switch control logics is 
shown in Table (3.2). Fig. (3.3(a)) illustrates the GMOS 
logical circuit for realizing this truth table. While four 
of the resistors are equal and of value R each, the fifth 
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3.1 - The Generalized Iraniccance Converter (GIC) 

Implementation Using OA' s 
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Filter 

Type 


h 


"2 


"3 
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(S^+W ^)/D(s) 
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c 
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QP 


Tl= is^- 


^ 2 

Qp s +wp ) /D (i 



where T(s) = N(s)/D(s) and D(s) = + (Wp/Qp)S + Wp^ 



TABLE 3.1. - The Elements Identification for Different 
Realizations of the GIC Filter . 
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Fig. 3.2 - Schematic Diagram of the Programmable GIC 
Filter Showing the Controlled Nodes 
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f'lg. 3.3 



Different Elements Realizations 
Corresponding Switches Used for 
Selecting the Filtering Type . 



and the 
Digital 
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TABLE 3 .2 -The Truth Table of the Switches Logic Used to 
Select the Filtering Function 
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30 . 



z Co 

^ c^- 




Fig. 3.3a The CMOS Logic Diagram Used to Control the 

Analog CMOS Switches of Fig. and to 

Realize the Truth Table 
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resistor is the Qp determining resistor and of value 
Rq=RQp. The two capacitors are equal and of value C=l/WpR. 
each. The two equal banks of capacitors are used to control 
ojp. Each bank contains n binary weighted capacitors 
connected in series through analog CMOS switches as shown in 
Fig. (3.5). Using a digital binary word of n bits to control 
Wp, 2^ different values of C will result at the 2 terminals 
of both capacitors banks that correspond to 2^ different 
values of p. Using a similar technique the value of Rq can 
be controlled through a bank of m binary weighted resistors 
in series, though analog CMOS switches as shown in Fig. (3.6). 
Using a digital binary word of m bits to control Qp, 2*^ 
different values of Rq can be achieved that correspond to 2^ 
different values of Qp. Thus, full independent control of 
the pole pair (bp and Qp are achieved by programming the 
switches to obtain the corresponding C and Rp. It can be 
easily shown that with minor modifications, an additional 
programmable element can be added for the control of the 
notch frequency. 

C. THE REALIZED GIC PROGRAMMABLE FILTER 

A complete circuit diagram of the constructed GIC filter 
is shown in Fig. (3.4). The values of m and n were selected 
to m=n=4. Thus, 15 different values of oJp (fp) and Qp were 
obtained as it is illustrated at the corresponding Table 
(3.3) and (3.4). the (designed) banks of the resistors for 
the control of Qp and the capacitors for the control of <-0p 
(fp) along with their control switches are shown 

correspondingly in Fig. (3.7) and Fig. (3.8). 
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Fig. 




3.4 - The Complete Circuit Diagram 

of the Programmable GIC Filter 
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digital word controlling Wp 



ccq Qj CX2 

0 0 



r 

'-'0 



1 


1 

i 


i 

> 

1 


1 








f 




^2- 


# 

( 

( 


1 

i 

i 


1 i 

j 


1 






-- 1 




V/ — 0/ 



Fig. 3.5 



The Two Capacitor Banks Realizations for the 
Programming of 
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The Resistor Bank Used to Realize R Needed 
for the Programming of Q^. • ^ 
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switch 
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0 


1 


6.34 


15.709 


0 


0 


1 


0 


7.16 


13.887 


0 


0 


1 


1 


7.72 


12.875 


0 


1 


0 


0 


8.42 


11.810 


0 


1 


0 


1 


9.20 


10.856 


0 


1 


1 


0 


10.00 


10.043 


0 


1 


1 


1 


11.00 


9.047 


1 


0 


0 


0 


15. 10 


6.650 


1 


0 


0 


1 


17.60 


5.637 


1 


0 


1 


0 


20.80 


4.823 


1 


0 


1 


1 


26.00 


3.828 


1 


1 


0 


0 


35.50 


2.805 


1 


1 


0 


1 


55.00 


1.810 


1 


1 


1 


0 


100.00 


.995 



Table 3.3 The four - fait words that control 
Fp and the corresponding capacitor 
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sw i tch 
control 


Rq 


Qp 


Sa 


Sb 


Sc 


Sd 


K, 


0 


0 


0 


0 


24.0 


15 


0 


0 


0 


1 


22.4 


14 


0 


0 


1 


0 


20.8 


13 


0 


0 


1 


1 


19.2 


12 


0 


1 


0 


0 


17.6 


1 1 


0 


1 


0 


1 


16.0 


10 


0 


1 


1 


0 


14.6 


9 


0 


1 


1 


1 


12.8 


8 


1 


0 


0 


0 


11.2 


7 


1 


0 


0 


1 


9.6 


6 


1 


0 


1 


0 


8.0 


5 


1 


0 


1 


1 


6.4 


4 


1 


1 


0 


0 


00 

• 


3 


1 


1 


1 


1 


3,2 


2 


1 


1 


1 


0 


1.6 


1 



Table 5.4 The four - bit words that control 
Rp and Qp 
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= lOOuf 
C 2 = 55uf 

= ■25a.F 

p = A ^ '-•-T 
^4 

Fig . 3.7 



The Constructed Blocker 
Capacitor Controlling the fp. 
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^ DIGITAL CONTROL] 



0 

0 


0 






t 

1 


T 

1 


r 

1 


f 

1 

1 


"1 


1 

..1 

1 


1 

1 

1 


1 

1^1 

1 

1 


1 

- 

1 — 


1 

J^oo 


I 

1 


1 

1 


-’AAA- 


AW 


AW 


AW- 


R, 


A 


^3 


A 



= 1.6k 
R 2 = 3.2k 
R^ = 6.4 k 
R^ = 12.8k 

Fig. 3.8 - The Constructed Block of 

Resistance Controlling Op. 
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■ IV . COMPUTER SIMULATION OF THE PROGRAMMABLE FILTER 



A. INTRODUCTION 

In order to observe the theoretical responses of the 
different realizations as well as to compare them with the 
experimental measurements, two computer programs were written 
in Fortran. those programs are shown in Appendices A and B. 
The first one simulates the programmable filter's responses 
Tl(s), T2(s) and T3(s) as functions of the network's 
admittances Yi, i=l . . .8, and the "constant" Op. Amps, gains 
Al and A2. A realistic value of A2=A2=10^ was given to the 
above gains corresponding to 741 Op. Amps, used in the 
experiment. the second program simulates the programmable 
filter's responses Tl(s), T2(s) and T3(s) as functions of the 
network's admittances Yi, i=l ... 8 and the frequency 
dependent Op. Amps, gains Al and A2. In this second case a 
single pole approximation value of wils was assigned for 
Ai(i=l, 2). Where, Wi is the GBWP of the Op. Amps, used, a 

g 

value of W1=W2=2“X xlO was given to the above Wi(i=l, 2) 
corresponding to 741 Op. Amps, used throughout the research. 

To be able to compare the computer simulation results 
with the experimental ones obtained, the same values of 
admittances were used as input data. That means that the 
same values of R, Rg and C according to Chapters III and V 
were used. With the value of R selected to be 1.6k Table 
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(4.1) illustrates tlie different values of Rqs used in the 
simulation to realize the different values of Qp and Table 

(4.2) illustrates the different values of capacitors used in 
the simulation to realize the different frequencies Wps. 

The "DO CASE I" command of Fortran simulated the digital 
logic (including the control switches) used to realize the 
different responses of the programmable filter that is LP, 
HP, BP, N and AP. The frequency's (wp) translations and the 
different QP values were simulated by changing the values of 
Rq and Cs at every run of the program. 

The above programs also simulated the transfer function 
of two cascaded GIC programmable filters as it will be 
discussed later in Chapter VI. Each of the filters could 
have been at different realization (as well as at different 
p and Qp) relative to the other one. The 25 possible 
combinations of transfer function realizations are shown in 
Table 4.3. 

B . SIMULATION RESPONSE ( S ) 

1 . Low Pass (LP) Realization 

Using the elements values prescribed in Table (3.1) 
yields the following L.P. transfer function; 



\2Cs-) 



”2 liip 

Q p 






(4.1) 
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For R=1.6k5 


Rq 

kS 


Qp 


1.1 


0.7 


1.6 


1.0 


3.2 


2.0 


6.4 


4.0 


11.2 


7.0 


16.0 


10.0 


32.0 


20.0 



Where Qp = ^ 



Table 4.1 



The Resulting Qp Values 
Introduced Rq Values in 
Simulation. 



for the Different 
the Computer 
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For R-1.6k2 



c 

nf 


fp 

khz 


■ 100.00 


.099 


50.00 


1.99 


26.00 


3.83 


12.5 


7.96 


6 . 6 


15.10 


3.3 


30.16 


0.85 


117.09 


0.42 


236.67 



fp = 1 

2 K.- R.C. 



Table 4.2 - The Resulting p Values for the Different 
Introduced C Value (R constant 1.6k?) 
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FIITEn 1 


filter 2 


J 


I 


• IP 


L? 


2 


2 


'HP 


HE 


1 


1 


(BP , 


BP 


4 


4 


N 


N 


3 


3 


- A F 


A? 


3 


3 


IP 


HP 


2 


1 


LP 


9P 


2 


4 


LP 


iN 


2 


3 


1? 


AP 


2 


3 


HP 


LP 


1 


2 


- HP 


BP 


1 


4 


HP 




1 


3 


H? 


A? 


1 


3 


BP 


L? 


M 


2 


BP 


HP 


4 


1 


3? 


N 


4 


3 


BP 


AP 


a 


•5 

J 


H 


IF 


3 


2 


N 


BP 


3 


1 


N 


BP 


3 


4 


H 


AP 


3 


3 


AP 


IP 


3 


2 


AP 


HP 


3 


1 


AP 


BP 


3 


4 


AP 


N 


3 


3 



ToWe. - The 25 Possible Combinations of 

Transfer Function Realizations 
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This transfer function assumes the following complex 



value 



>2(1^430- -QjQp 



(4.2) 



with magnitude of 






(4.3) 



Fig. 4.1 illustrates the theoretical "ideal" LPF 
magnitude response for fp=3.8Khz and for 3 different values 
of Qp. The simulation results match the equations of (4.2). 

Fig. 4.2 illustrates the theoretical "ideal" LPF 
magnitude response for Qp=2 and for 3 different values of^'p. 
the Op. Amps, gains Ai, (i=l, 2) are frequency depended. 
This dependence affects the magnitude of the filter and 
causes a frequency shift from the theoretical value of the 
ideal's case. Figs. (4.3) and (4.4) illustrates the ideal vs. 
nonideal theoretical LPF amplitude responses. 

Data extracted from Figs. (4.3) and (4.4) are 
illustrated in Tables (4.4(a)) and (4.4(b)) simultaneously. 
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AMPUTUDE D6 



LP.F AMPLRESPONSE (P=3.828KH) 




Fig. 4.1 - "Ideal" LPF Amplitude Response 
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AMPLITUDE DB 



O 



LP.F AMPLRESPONSE (Q=2) 




Fig. 4.2 - "Ideal" LPF Amplitude Response 
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AMPLITUDE DB 



LP.F AMPLRESPONSE (Q=2) 




ideal 

nonideal 

Fig. 4.3 Ideal .vs. nonideal L.P.F. amplitude response 
and frequencies (1.99K, 7.96K, 30.16K and 117 



for Q=2 
lK)h> 
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AMPUIIUDE DB 



LP.F AMPUTUDE RESPONSE 




2. Q = 2.0 

3. Q = 4.0 

4. Q = 7.0 

Fig. 4.4 Ideal .vs. nonideal L.P.F. amplitude responce for 
frequencies (7.96K, 30.lK)Ktand variety of Qs. 
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Data From Fig. (4.4) LPF Characteristic 





Ideal 




Nonideal 


QP 




1T2(s)| 
picJc value 


Rq 


QP 


|T2(s) ! 




at 






at picK value 


.07 




3.0 


1.1 


6.7 


3.0 


2.0 




12.0 


3.2 


2.0 


12.0 1 


4.0 




18.1 


6.4 


4.0 


18.1 I 


7.0 




22.9 


17.2 


7.0 


22.7 


The piclc 


values obtained 




The pic)c values obtained 


at 


fpiclc = 7.9)c 




at fpiclc - 7.9)c 

1 


(b) 




Ideal 


C 


? 

1 

Nonideal 


' Fp 




Amp. Response 


nf 


fp 


Amp. Response 


)ch2 




dB 




]<h2 


: dB 

1 


1990 


.45 


12.04 


50 


1990.45 12.72 


7165 


.60 


12.24 


12.5 


8359.87 12.113 


28662 


.42 


12.31 


6.6 


25875.80 12.04 


109745 


.50 


12.31 


0.085 


90565.00 11.08 



Table 4.4 - Data from Figs. (4.3) and (4.4) indicating 
the Affect of Frequency Dependency of 
Ai , i=l , 2 
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2 . High Pass (HP) Realization 



Using the elements value shown in Table (3.1) yields 
the following HP transfer function; 






Qp 



takes the following complex value at wp. 



which magnitude is 

I '^GoDp) 1 - '2G}pz ^ <3b 



(4.4) 



(4.5) 



(4.6) 



Fig. (4.5) illustrates the theoretical "ideal" HPF 

amplitude response for fp=3.8KHZ and for 3 different values 
of Qp. The computer simulation results match that of the 
(4.6) relation. 

Fig. (4.6) illustrates the theoretical "ideal" HPF 

amplitude response for Q=2 and for 3 different frequencies. 

Fig. (4.7) illustrates how this dependency effects 

the amplitude of the HP filter. (The amplitude decreases as 
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AMPLITUDE DB 



H.P.F AMPL RESP0NSE(P=3.828KH) 




Fig. 4.5 - "Ideal" HPF Amplitude Response 
for a Variation of Qs . 
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AKIPUTUDE DB 



H.P.F AMPLITUDE RESP0NSE(Q=2) 




Fig. 4.6 - "Ideal" HPF Amplitude Response. 
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AMPUTUDE DB 



HJ>.F AMPUTUDE RESP0NSE(Q=2) 




Fig. 4.7 -Nonideal HPF Amplitude Response. 
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the frequency increases) while Figs. (4.8) and (4.9) 

illustrate the frequency shift of the amplitude response due 
to that dependency. 

From the date of Fig. (4.8), Table (4.5) were 

constructed. 

3 . Band Pass (BP) Realization 

Using element values from Table 3.1 the following BP 
transfer function is achieved; 




(4.7) 






talces the following value atujp. 






(4.8) 



which has constant magnitude of 6dB 




(4.9) 
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AMPLITUDE DB 



H.P.F AMPLITUDE RESP0NSE(Q=2) 




Ideal 

nonideal 



Fig. 4.8 - "Ideal" vs. Nonideal HPF Amplitude 
Response . 
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Ideal 



c 



Nonideal 



fp Amp Response 

Hz dB 


nf 


fp 

Hz 


Amp Response 
Db 


4,378.98 


12.28 


26 


4,378.98 


12.148 


8,359.87 


12.30 


12.5 


8,359.87 j 


i 12.11 


15,923.57 


12.31 


6 . 6 


15,923.57 ' 


' 12.02 


31,847.13 


12.31 


3.3 


31,847.73 j 


1 11.41 


123,407.00 


12.31 


0.85 


107,484.00 j 


1 9.21 


252,781.05 


12.31 


0.43 


187,101.80 1 


' 9.21 1 


242,781.05 


12.32 


0.43 


187,101.80 


6.69 1 

1 



Table 4.5 - Data Illustrating the Ideal vs. the Nonideal 
Responses of the HPF 
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AMPUTUDE DB 



H.P.F AMPL RESP0NSE(P=8K/16K) 




2. Q = 2.0 

3. Q = 4.0 

4. Q = 7.0 

Fig. 4.9 - "Ideal" vs. Nonideal HPF Amplitude 
Response . 
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Fig. (4.10) illustrates the theoretical BPF magnitude 
response for fp=3.83 KHZ for different values of Qp. These 
agree with the (4.9) equation since as it is indicated by the 
simulation plot the amplitude is constant and independent of 
Qp. 

Fig. (4.11) illustrates the above concept but at 
fp=15.1KHZ. 

Fig. (4.12) and (4.13) illustrate the theoretical 
"ideal" BPF amplitude response for Qp=2 and for different 
values of p. As it is indicated from Fig. (4.13) the 

amplitude remain constant even at very high frequencies (10 
HZ). But with Ai, (i=l, 2) depending on frequency the 
amplitude decreases as the frequency increases. This is 
indicated in Figs. (4.14) and (4.15) which describe the BPF 
amplitude response plots for Q=4 and Q=1 respectively and for 
different frequencies. The frequency dependence of Ai, (i=l, 
2) creates a frequency shift from the ideal theoretical value 
which is indicated in Figs. (4.16) and (4.17). Table (4.6) 
illustrates the data extracted from figs. (4.16) and (4.17). 

4 . Notch (N) Realization 

Using the admittances value of Table (3.1), the 
following Notch transfer function can be achieved: 



(4.10) 
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AMPLITUDE DB 



B.P.F AMPL RESP0NSE(P=3.828KH) 




FREQUENCY HZ 



Fig. 4.10 - "Ideal" BPF Amplitude Response. 
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AMPLITUDE DB 



B.P.F AMPLITUDE RESP0NSE(Fi=16K) 




Fig. 4.11 -"Ideal" BPF Amplitude Response. 



90 



AMPLITUDE DB 



B.P.F AMPUTUDE RESP0NSE(Q=4) 




Fig. 4.12 - ’Ideal" BPF Amplitude Response. 
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AMPLITUDE DB 



B.P.F AMPUTUDE RESP0NSE(Q=4) 




Fig. 4.13 - "Idea]!' BPF Amplitude Response. 



AMPUTUDE DB 



B.P.F AMPLITUDE RESP0NSE(Q=4) 



O 



o 

n' 



q 



j 




I '’ f 7 T l “ l r 

irf 






T 1 I I r 



n 



REQl^NCV 



Fig. 4.14 - Nonideal PBF Amplitude Response. 
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AMPUTUDE DB 



B.P.F AMPUTUDE RESP0NSE(Q=1) 




Fig. 4.15 - Nonideal BPF Amplitude Response. 
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AMPUTUDE DB 



B.P.F AMPUTUDE RESP0NSE(Q=4) 




Ideal 
non ideal 



Fig. 4.16 - Ideal vs. Nonideal BPF Amplitude 
Response . 
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AMPUTUDE DB 



B.P.F AMPLITUDE RESP0NSE(P=16K) 




Ideal 

nonideal 



Fig. 4.17 - Ideal vs. Nonideal BPF Amplitude 
Response . 
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Ideal 


C 

nf 


Nonideal 


fp Amp. 

khz 


Response 

dB 


fp Amp. 

khz 


Response 

dB 


7,961.78 


6.02 1 


12.5 


7,969.738 


5.958 


30,254.78 


6.02 ] 


3.3 


28,662.42 


5.75 


115,445.80 


5.97 ! 

1 


0.85 

1 


95,541.38 


4.92 



and 



Ideal 


! Rq 

k 


Nonideal 


QP 


Max. 

Amp. Response 
dB 


QP 


Max. 

Amp. Response 
dB 


4 


6.01 


6.4 


4 


5.902 


8 


6.01 


12.8 


8 


5.728 


20 


5.95 


32.0 


20 


5.528 


fp^ = 15, 127. 33Hz 




fPl 


= 14, 739. 29Hz 



Table 4.6 - Data Indicating the Ideal vs. the Nonideal 
Responses of the BPF Realization 
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AMPLITUDE DB 

38.0 -« 7.0 - 81.0 -160 - 0.0 - 3.0 



N.F AMPUTUDE RESPONSE(Q=2.) 




Fig. 4.18 



Ideal Notch Amplitude 



Response . 
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NOTCH AMPLRESP0NSE(P=3.83KH) 




Fig. 4.19 - Ideal Notch Amplitude Response. 
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AMPLITUDE DB 



N.F AMPUTUDE RESPONSE(Q=S.) 




Ideal 

nonideal 



Fig. 4.20 - Ideal vs. Nonideal Amplitude Response. 



lOG 



where n is the Notch frequency. 



At P the transfer function takes the value 






G?p 




(4.11) 



Fig. (4.18) illustrates the "ideal" theoretical Notch 



filter amplitude response for a variety of frequencies and 
constant Qp(Q=2), while Fig. (4.19) for a variety of Qps and 
for constant p ( fp=3 . 83KHZ) . 

Fig. (4.20) illustrates the effect of the frequency 
dependency of Ai, (i=l, 2) on the response. 

5 . All Pass (AP) Realization 

As proposed in Table (3.1) using the same admittances 
values as Notch. An All Pass transfer function can be 
derived as 




(4.12) 




qp 
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which takes the following values 
o\j iD 0 , \Ti ^ 4 ~0 y 

UO-^oo^ 1 ~i(joo)\ — ^2 = 6ab^ 



/ >t Cjo") ^ 

LZlSi^ -> o° 



(4.13) 



The above agree with the computer simulation results 
of Figs. (4.21) and (4.22). 
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AMPUTUDE DB 



AJ> AMPUTUDE RESP0NSE(Q=2) 




Fig. 4.21 - Ideal vs. Nonideal APF Amplitude 
Response . 
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PHASE DEGREE 

e0.0 ISO.O 180.0 840.0 300.0 380.0 



A.P PHASE RESPONSE(Q=2/P=8K) 




Fig. 4.22 - Ideal vs. Nonideal Phase Response. 
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V. REALIZATION OF PROGRAMMABLE GIC FILTER 



A. EXPERIMENTAL RESULTS 

After the circuit of Figs. (3.10), (3.9), (3.8), (3.7), 

was constructed a variety of measurements were ta)<en in order 
to study the response of the networlc to the different inputs 
(control bitwords). To observe the affect of the control 
switches which introduce a resistance of 80 each at CLOSED 
position, two values of Rs, Rq, and Cs were used with one 
decade difference in magnitude. That means that R was given 
the values of 1.6K (as discussed in Chapter III) and 16K, the 
four resistors that consisted the Rq ban)c were of values 

(1.6K, 3.2K, 6.4K, 12. 8K) in the first case and (16K, 32K, 

64K and 128K) in the second one, and that the capacitor 

banlc's capacitors were chosen of values (lOOnF, 50nF, 12, 

5nF, llnF) and (lOnF, 5nF, 1.2nF, O.lnF) accordingly, to be 
able to )ceep the range of frequencies as much the same as 
possible for both cases. 

1 . Low Pass Filter 

With the topology-control bit word 000 the networlc 
realized a LPF response. Fig. 5.1(a) illustrates the LP 
response for R=1.6Kito a variety of frequencies for Q=5 
while 5.1(b) illustrates the same but for Q=2. It can be 

observed from 5.1(a) (Q=5) that the amplitude response 
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(a) 




1 - f = 0.99khz (b) 

2 - f = 1.98khz 

3 - f = 3.83khz 

4 - f = 7.96khz 

5 - f = 12.0khz 



Fig . 5.1 



(a) LPF Response for R=1.6k (Q=5, Rq=8.0k) 

(b) LPF Response for R=1.6k (Q=2, Rq=3.2k) 
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decreases while the frequency increases as it was expected 
from computer simulation. This is occurred up to the 
frequency of 6KHZ; then it started increasing with the 
frequency, while at Fig. (5.1(b)) (Q=2) it remained constant. 
Fig. (5.2) illustrates (for R=16KHZ) that the observation in 
Fig (5.1) is not any more the case (for that frequency range) 
and the networ)< responses the same as in computer simulation 
while in Fig. (5.3) (Q=5) the above can be noticed again, 

this is due to the interference of the control switches as it 
was discussed previously. Figs. (5.4) and (5.5) illustrate a 
variation of Q values for f=9KHZ for the two cases (R=1.6k 
and R=16K ) respectively. A difference in magnitude can be 
observed due to the interference of the control switches. 
Figs. (5.6) and (5.7) illustrate the same but for F=12.8KHZ. 
Figs. (5.8) and (5.9) illustrate the phase and amplitude 
response of the LPf. 

2 . High Pass (HP) Realization 

With the topology-control bitword 001, the network 
realized a high pass filter. Figs. (5.10) and (5.11) 

illustrate the HPF amplitude response for a variety of 
frequencies and (Q=2). It can be observed (as in LPF 
realization) that for this frequency range and for R=1.6Kt 
the amplitude starts increasing as the frequency increases, 
decreases then again while for R=16KT this does not occur, 
the same observation submerges comparing Fig. (5.14) and 
(5.15). Figs. (5.12) and (5.13) are the plots obtained for 
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1 - f = 0.99khz 

2 - f = 1.99khz 

3 - f = 3.83khz 

4 - f = 7.90khz 

5 - f = 12.8khz 

Fig. 5.2 - LPF Amplitude Response 
R=1.6k (Q=2, Rq=32k) 
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2 - f = 1.99khz 

3 - f = 9.20khz 

4 - f = 7.96khz 

5 - f = IS.Okhz 



Fig. 5.3 - LPF Amplitude Response 
R=16k (Q=5, Rq=50k) 
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1 - Q = 1 

2 - Q = 2 

3 - Q = 4 

4 - Q = 7 



Fig. 5.4 - LPF Amplitude Response (R=1.6k), 
f=7.96khz 



110 




1 - Q = 1 

2 - Q = 2 

3 - Q = 4 

4 - Q = 7 



Fig. 5.5 - LPF Amplitude Response (R=10k) , 
f=7.96khz 
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2 - Q = 2 

3 - Q = 4 

4 - Q = 7 

Fig. 5.6 LPF Amplitude Response (R=1.6k) and 
f=12.8khz 
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1 - Q = 1 

2 - Q = 2 

3 - Q = 4 

4 - Q = 7 

Fig. 5.7 - LPF Amplitude Response (R=1.6k), 
f=12 . 8khz 
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1. Phase Response. 

2. Amplitude Response. 

Fig. 5.8 - LPF Amplitude/Phase Response. 
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Fig. 5.9 - LPI Phase/Amplitude Response. 
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1 - f = 0.99khz 

2 - f = 1.99khz 

3 - f = 3.83khz 

4 - f = 7.96khz 

Fig. 5.10 - HPF Amplitude Response 
(R=1.6k), Q=2 
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1 - f = 0.99khz 

2 - f = 1.9 khz 

3 - f = 5.0 khz 

4 - f = 8 . 0 khz 

5 - f = 13.8 khz 



Fig. 5.11 - HPF Amplitude Response 
(R=1.6k), Q=2 
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1 - Q = 0.7 

2 - Q = 2.0 

3 - Q = 4.0 

4 - Q = 7.0 



Fig. 5.12 - HPF Response (R=1.6k), 
f=6 . 65khz 
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1 - Q = 2.0 

2 - Q = 4.0 

3 - Q = 7.0 



Fig. 5.13 - HPF Amplitude Response 
(R=1.6k), f=6.65khz 
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1 - f = 0.99khz 

2 - f = 1.99khz 

3 - f = 3.83khz 

4 - f = 7.96khz 

Fig. 5.14 - HPB Amplitude Response 
(R=l . 6k) , Q=5 
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1 - f = 0.99khz 

2 - f = 1.99khz 

3 - f = 3.83khz 

4 - f = 7\ 96khz 

5 - f = 12.8khz 



Fig. 5.15 - HPB Amplitude Response 
(R=1.6k), Q=5 
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f=6.65KHZ and a variety of Qs for R=1.6K4 and R=16K2 , 



respectively. a difference of about 12dB in amplitude 
appears. Finally, figs (5.16), (5.17), (5.18), (5.19) 
illustrate the phase and amplitude responses of the HP 
realization. 

3 . Band Pass (BP) Realization 

The 010 topology-control bitword realizes the Band 
Pass Filter. Figs. (5.20) and (5.21) illustrates the 
amplitude response for set of frequencies and Q=10. At both 
cases the amplitude decreases while the frequency increases 
until the value of 9KHZ. Then starts increasing again. It 
can be also observed that for R=1.6K the frequency deviation 
from the theoretical fp is larger. 

Figs. (5.22) and (5.23) illustrate the amplitude 
response again, but for Q=l. This time the deviations from 
the theoretical response (both of amplitude fluctuations and 
frequency shift) are less. 

Figs. (5.24), (5.25), (5.26), and (5.27) are also 
plots of the amplitude response but for variation of Q. A 
difference of approximately IdB appears for the lowest 
frequency (3.8KHZ) and of 0.5dB for the higher (9KHZ), 
Finally, Figs. ( 5 . 28 ) - ( 5 . 32 ) illustrate the phase and 
amplitude of the BPF response. 

4 . Notch (N) Filter Realization 

With the topology-control bitword Oil, a Notch filter 
realization can be achieved. Figs. (5.33) and (5.34) 



122 



XFR FC 
X F R F C 



H M K R : 2 . 7 d B 

H ' 0 CENTER 



2 dB/Diy 
5 0 D I 



s 


















' 


\1 


- 




















J. 






2 






1 — 

< 

r* 

l_ 
































































































1 

Li 























•'. 0 Hz 

M K R ' 1 4 0 0 H : 



2 5 f •: H z . 
BH = 20 0 H: 



ux) 



X F R F C T N M K R > 2.9 d B 2 d B ■ - D I 

X F R F C T H : 0 CENTER 5do . - D I \> 























f 


s, 

1 u 


i‘ 
















I 


/ 

1 




















l' 

f 






- — — 















rT 




















/ 


















L_J 


/ 


















r 



















\ 0 H z 2 5 K H z / 

M K R : .3 800 Hz B W = 2 0 0 H z 






1. Phase 

2. Amplitude 

Fig. 5.16 - HPF Amplitude/Phase Response 
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Fig. 5.17 - HBF Amplitude/Phase Response 
(Q=0.5) 



124 



XFR FC 
X F R F C 



U M K R : 3 . 6 d B 

M = 0 0 CENTER 



2 d B D I W 
5 0 o / D I 





___ 


r'-- 

/ 




^Ul 










! 


















r. .V . 




1 


















V 














__1 


___ 








1 




















\ 

1 






































L . 


1 



















0 Hz 

MKR ^ A 800 Hz 



25 KHz , 
B W •• 2 0 0 H : 



X F R F C 
X F R F C 



N M K R •• 2 b d B 2 d B • D I V 

N • 0 o C E N T E R 5 0 o D I V 























\ 


i 


V“» 






2. 














\ 






r- 
















1 























































































il 























0 Hz 

M K R = 1 2 0 0 Hz 



2 5 K H z / 
B W : 2 0 0 Hz 



1. Phase 

2. Amplitude 



Fig. 5.18 - HPF Phase/Amplitude Response 
( 0 = 2 ) 



125 



X F R F C T H t1 K R •• 3 . 4 d B 2 d B D I 'J 

XFR FCTN= 0O CENTER 50O/DIV 





V { 


V 


















i 




■H- 


1 . ^ 












1 


1 


V 


















/ 






. j. 





... 






























1 




















1 


















1 / 




















\ 0 
n K R ^ 


Hz_ 


5 600 H 






2 5 
B N = 


1 KH 

2 00 


T 

“Hz 



1. Phase 

2. Amplitude 



Fig. 5.19 - HPF Amplitude Response (Q=2) 



✓ 



126 




1 - f = 0.99khz 4 - f = 6.65khz 

2 - f = 1.99khz 5 - f = 10.5khz 

3 - f = 3.8khz 6 - f = IS.lkhz 



Fig. 5.20 - BPF Amplitude Response 
(R=1.6k), Q=10 
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1 . 

2 . 

3. 



Fig. 5,21 - BPF Amplitude Response 
(R=1.6k), Q=10 



f = 0.99K 
f = 1.99K 
f = 3.8K 



4. f = 7.96K 

5. f =11.2 K 

6. f =15.1 K 
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1. 


0 . 99khz 


4. 


6 . 65khz 
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1.99khz 


5. 
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6. 


15 . Ikhz 



Fig. 5.22 - BPF Amplitude Response 
(R=1.6k) , Q=7 
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1. 0.99khz 

2. 1.99khz 
3 . 3 . 83khz 



4. 7.96khz 

5. 12.8khz 

6. 15. khz 



Fig. 5.23 - BPF Amplitude Response 
(R=1.6k), Q=1 



130 




2 - Q = 4.0 

3 - Q = 1.0 



Fig. 5.24 - BPI Amplitude. Response 
(R=1.6k), f=3.83khz 
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2 - Q = 4.0 

3 - Q = 1.0 



Fig. 5.25 - BPI Amplitude Variation 
(R=1.6k), f=3.83khz 
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1 - Q = 10.0 

2 - Q = 4.0 

3 - Q = 1.0 



Fig. 5.26 - BPF Amplitude Response 
(R=1.6k) 
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1 - Q = 10 

2 - Q = 4 

3 - Q = 1 



Fig. 5.27 - BPF Amplitude Response 
(R=l . 6k) , f=9 .khz 
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Fig. 5.33 - Notch Amplitude Response 
(R=1.6k), Q=4 
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1 - f = 0.99khz 

2 - f = 1.99khz 

3 - f = 3.83khz 



Fig. 5.34 - Notch Amplitude Response 
(R=1.6k) 
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illustrate the amplitude response for different frequencies 
and for Q=4, while Figs. (5.35) and (5.36) illustrate the 
amplitude response for f=lKHZ and a variety of Qs. Fig. 
(5.37) illustrate the phase response in addition to the 
amplitude one. 

5 . All Pass (AP) Realization 

The topology-control bitword 100 realizes the 
All-Pass filter. Figs (5.38), (5.39) and (5.40) illustrate 
the plotted amplitude and phase responses. 

B . CONCLUSION 

The constructed circuit performed as predicted by the 
theoretical analysis and the computer simulations. This 
means that it realized all the desired filtering transfer 
functions LP, HP, BP, N, and AP. The effect of interference 
of the control switches nonideal performances which is more 
severe at high Qps can be minimized by increasing the values 
of R (meaning that the values of Rq also increases and the 
values of Cs decreases). 
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3 - Q = 8.0 



Fig. 5.36 - Notch Amplitude Response 
(R=1.6k) , f=0.99khz 
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2 - Q = 4.0 

3 - Q = 8.0 



• Fig. 5.37 - Notch Amplitude Response 
(R-1.6k), f=0.99khz 
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2 - Q = 8.0 

3 - Q = 12.0 



Fig. 5.38 - APF Amplitude Response 
(R=1.6k), f=2.8khz 



146 




Fig. 5.39 - APF Amplitude Response 
(R=1.6k) for f=0.99khz, 
f=1.99khz, and f=2.8khz 
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VI. COMBINING HIGHER ORDER SECTIONS 



By cascading two or more programmable filters, higher 
order transfer functions can be obtained. Fig. (6.1) 
illustrates two cascaded GIC programmable filters. 




functions Ti(s) and Tj(s) can obtained as it is indicated in 
Table (4.3). 

For Lp-LP combination and for ideal theoretical case 
(Al=A2-->oo) a fourth order low pass filter can be obtained 
with transfer function: 

2 






( 6 . 1 ) 



The computer simulation of the fourth order transfer 
function of a nonideal theoretical filter vs. the second 
order case is illustrated at Fig. (6.2). The experimental 
results taken from the properly designed and built circuit as 
illustrated in Fig. (6.3) are shown in Fig. (6.4). 

Using the same procedure as above a fourth order HP-HP 
combination for the ideal theoretical response is given by 



\aC^').= 4 





S V 




which at Wp takes the complex value of 



( 6 . 2 ) 



(6.3) 



with magnitude of 



I \ i 1 = A 0^ = 4 O i og t 



(6.4) 



Fig. (6.5) illustrates the fourth order nonideal 
theoretical HP filter response vs. the second order one. 
Both at Qp=2 and fp=8KHZ , while fig. (6.6) illustrates the 
experimental responses . For Qp=2 according to (6.4) and 
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Fig. 6.2 - LPF Fourth Order vs. Second Order 

Ideal Response from Computer Simulation. 
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Fig. 6.3 - Network with Logic for Realizing 
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AP Transfer Functions. 
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1. single 

2 . Cascade 



Fig. 6.4 - 4th Order vs. 2nd Order 

Experimental LPF Response 
(Q=2) 



153 



AMPLITUDE DB 



H.P.F AMPLRESP0NSE(Q=2/P=4K) 




Fig. 6.5 - "Ideal" 4th vs. 2nd Order HPS Responses 
From Computer Simulation. 
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Fig . 6.6 



4th Order vs. 2nd Order HPF obtained 
Responses from the Constructed Circuit. 
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(2.4) a difference of 12dB were expected between the second 
and fourth order filter which is approximately the case in 
both computer simulation and experimental results. 

The BP-Bp combination leads to the theoretical "ideal" 
fourth order transfer function: 




which takes the following value at u)p: 



(6.5) 



( 6 . 6 ) 

According to this, a magnitude response of 12 dB 
approximately at tCp, was expected from both experimental and 
computer simulation results. Fig. (6.7) illustrates the 
simulated response of a fourth order nonideal HP filter vs. a 
second-order one. A difference of 5.8 dB instead of 6 dB can 
be observed. Fig. (6.8) illustrates the experimental 
response, where a difference of 5.4 dB appears basically due 
to the interference of the control switches and the 
nonideally matched values of the capacitors which control the 
fp selection. 
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AMPUTUDE DB 
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— , — . — cascade 
single 



Fig. 6.7 - Fourth Order vs. Second Order BPF 
Amplitude Responses from Computer 
Simulation . 



157 




Fig. 6.8 - Fourth vs. Second Order BPF 
Responses Obtained from the 
Constructed Circuit. 
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The N-N combination for the ideal theoretical case 



(Al=A2-->oo) results the following transfer function; 



T2C5')- ( 



(6.7) 




where ti>n is the Notch frequency. 

Fig. (6.9) illustrates the fourth order "nonideal" Notch 
filter simulated response vs. the second order one. 

The AP-Ap combination for the theoretical ideal case 
result in the following transfer function; 




( 6 . 8 ) 



which takes the following values for S=0 and oo; 



Tl i 



(6.9) 



with amplitude and phase of 






( 6 . 10 ) 
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Fig. 6.9 - Fourth Order vs. Second Order Notch 

Filter Amplitude Response from Computer 
Simulation . 
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and 



TiLjoo^ ^ 115 



( 6 . 11 ) 



with amplitude and phase of 

lTi(jcD^i_- -- (f ab 



The computer simulation of nonideal fourth order all pass 
filter vs. a second order one which is illustrating at Fig. 
(6.10) matches the above. 

Fig. (6.11) illustrates the resulting Chebychev filter 
from a BF-BP combination with Qp=4 and different frequencies 
(fpl=8KHZ fp2=10KHZ), while Figs. (6.12) and (6.13) 
illustrate the resulting Chebychev filters from the designed 
and built circuit. Fig. (6.14) illustrates the resulting 
response from a HP-LP combination. 
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AMPLITUDE DB 
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Fig. 6.10 - APF Fourth Order vs. Second Order 
Amplitude Response from Computer 
Simulation . 
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Fig. 6.11 - Chevycheb Response for f^=8k, 

f2=10k 



163 




Fig. 6.12 - Chevychev Response Obtained by 
Cascading BP-BP 
fj^=6.65khz £ 2 = 12 . 8khz 

Q^=3 
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Fig. 6.13 - BP-BP Response Qj^=Q2=4 
(f^=3khz, f2=15khz) 
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Fig. 6.14 - HP-LP Response Qj^=Q2=4 
(f^=3khz, 1^=9. Okhz) 
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VII. APPLICATION OF THE PROPOSED GIC PROGRAMMABLE FILTER 
IN FREQUENCY HOPPING SYSTEMS 

A. BACKGROUND 

1 . General Description of Frequency Hopping Signals 

Frequency hopping is a spread spectrum modulation 
technique used to generate many possible carrier frequencies 
over a large bandwidth. Of all the possible carrier 
frequencies, only one is selected at a given time. However, 
all frequencies are eventually selected during some time 
interval. 

Frequency hopping (FH) may be pictured as an RF 
carrier whose center frequency is "hopped" over many 
frequencies. The hopping may be either in a simple sequence 
or a pseudorandom sequence. 

The hopping rate of a frequency hopping system does 
not affect the bandwidth of the output spectrum. In a direct 
sequence system the chip rate determines the total bandwidth ► 
In a frequency hopping system, however, the bandwidth is 
determined by the highest and lowest frequencies of the 
frequency hopped carriers. For example, if the highest 
frequency carrier is at 15 MHz and the lowest frequency 
carrier is at 10 MHz, the total signal bandwidth is 5 MHz. 
This is the bandwidth regardless of the hopping rate. This 
allows wideband spread spectrum signal generation at low 
hopping rates. 
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2 . Signal Generation 

Frequency-hopped signals may be generated in several 
ways. The different methods are classified into two groups: 

(1) Direct synthesis, and 

(2) Indirect synthesis. 

One important aspect of frequency hopping synthesis 
is coherency. coherent signal synthesis is defined as the 
establishment of a known and repeatable phase each time a new 
frequency is hopped to. Non-coherent signal synthesis is 
defined as the establishment of a random or unknown phase 
each time a new frequency is hopped to. some techniques, 
direct or indirect, of signal generation can be used as a 
coherent frequency source. In other techniques, the changing 
of frequencies creates non-coherent sources. 

If a frequency hop system is a coherent, it will have 
a signal-to-noise advantage over a non-coherent system. 

a. Direct Synthesis 

The direct approach to signal synthesis utilizes 
techniques which enable direct synthesis of different 
frequencies. Examples of direct synthesis techniques are: 

(1) Frequency mixing, and 

(2) Surface acoustic wave devices. 

Frequency mixing for single synthesis is a common 
technique used to generate many different frequencies. An 
example of the frequency mixing technique is show in Fig. 
(7.2). 
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figure; 7.1 BLOCK DIAGRAM OF A FREQUENCY HOP MODEM 

A pseudorandom code generator selects one of many transmit frequencies during a small 
time interval. Tlie traffic modulates the carrier frequency which is spread over many 
different frequencies by the hopping action. The receiver dehops the input signal 
into a narrowband IF. The code synchronizer locks the pseudorandom code generator in 
the receiver to tfie received signal. A data demodulator removes the traffic from the 
IF anipl i f ier output . 
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In Fig. (7.2), an RF switch selects one of 
several frequency inputs. Two of these input signals of 
different frequencies are multiplied together to generate a 
new output frequency. The device used to multiply the two 
signals together is called a frequency mixer. when two 
frequencies are mixed, the sum and the difference of the 
frequencies are generated. In order to select only one of 
these frequencies, a "filter is used to reject the unwanted 
frequency." A filter tuned to the desired frequency would 
allow selection of that frequency while rejecting the other. 
By selecting the mixing frequencies in the proper order, the 
output frequency can be stepped through several different 
frequencies. At each frequency mixer output, a filters is 
required to reject unwanted frequencies. The filters may 
require a short time for the signal to stabilize after it is 
selected. The time required for the filter to stabilize at 
each new frequency may ultimately determine the maximum 
hopping rate of the direct frequency synthesizer, 
b. Indirect Synthesis 

The indirect method of signal synthesis is 
defined as frequency synthesis through .the use of 
phase-locked oscillators. One common indirect method for 
synthesis is shown in Fig. (7.3). 

In this circuit a phase-locked loop is used to 
generate the numerous carrier frequencies. The phase-locked 
loop has an internal oscillator whose output frequency, Fo, 
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riGURE 7.2 DIAGRAN OF A DIRECT FREQUENCY SYNTHESIZER 

A psciKiorarulum sequpnce <;clpct‘; diffcrpnt combinations of F| tijrough , which are mixed 
together to form a new frentiency. Each new frequency is mixed with the rf oscIMator 
for tran'mission as one of the many hopping frequencies. 
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FREQUENCY 

REFERENCE 




FIGURE 7.3 INDIRECT SIGNAL SYNTHESIS USING A PHASE-LOCKED LOOP 

When the reference frequency, F^ Is at the same frequency as F^*, the voltage- 
controlled oscillator produces a constant output frequency, Fq. Fq Is 

divided by _the programmable divider, Fq Is equal to N times F^ . By changing 
the divider ratio, N, many output frequencies are possible. 
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is shown in Figure (7.3) as the phase-locked loop output. The 
divide-by-N circuit divides this oscillator frequency by a 
selected number, N. 

The phase-locked loop internally adjust its 
output frequency Fo so that FI* is the same frequency as the 
reference frequency FI. If the divide-by-N circuit output 
frequency is initially lower than the reference frequency FI, 
the oscillator output frequency is automatically increased 
until FI and FI* are identical. When this occurs, the 
oscillator output frequency will become stable and remain at 
that frequency until the number, N, changes. When this 
number changes, the oscillator frequency is again 
automatically adjusted so that FI and FI* are again 
identical. 

B. PROPOSED USES OF PROGRAMMABLE FILTER 

A wide field of applications exist in FH systems for the 
programmable GIC filter. The outstanding performance of the 
filter (including sensitivity and stability) and its high 
speed response to the different inputs (due to the use of 
CMOS integrated circuits) make it very exceptional in this 
field. The first proposed use is indicated, in Fig. (7.4). 
The figure illustrates a receiver of a frequency hopping 
demodulator. The received frequency hopped signal after 
heterodyned by the RF mixer passes through the GIC filter in 
a BP realization (at this application the topology control 
network does not need to exist since the BP realization is 
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A Frequency 
Programmable 



Fig 



7.4 



• OSCI LLATOR Froc^uency 

Hop Clock 



the only type to be used) . The frequency shift of the filter 
is controlled by the synchronized pseudorandom code 
generator. This code generator also controls Q through some 
interfaced binary logic in order to correct the amplitude 
reduction which appears at high frequencies. The frequency 
shift problems can be easily corrected by the use of 
composate Op. Amp. Fig. (7.5) as it is extensively analyzed 
in Refs. [32], [34], and [35]. 

The programmable filter can also be used in the direct 
frequency synthesizer as it is illustrated at Fig. (7.2). At 
each frequency mixer output exists the need of a filter 
required to reject the unwanted frequencies. The frequencies 
to pass are not always the same, but they hop. the BP 
realization of the filter is used which center frequency can 
be controlled accordingly. The filter may require a short 
time for the signal to stabilize after it is selected. The 
time required for the filter to be stabilized at each new 
frequency may ultimately determine the maximum hopping rate 
of the direct frequency synthesizer. 
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Practical BP FiUer Realization of the 
posite GIC 



176 



VIII. CONCLUSION 



The novel design described here has resulted in a 
universal prograinmable filter than can be digitally 
controlled to realize almost any practical filter 
specifications. This is done through the use of CMOS 
switches controlled by binary codes to program the order of 
the filter, the filter topology, the filter center frequency 
and selectivity. The design procedure required developing 
optimum switching arrangements for the minimum redundancy in 
components and least dependence of the filtering function on 
switching imperfections such as switches stray capacitances 
and non-zero and nonlinear switch-on resistance. The 
sensitivities of Qp, Wp are found to be low with respect to 
the passive and active elements variations. The experimental 
result show close agreement between theory and practice. 
Further, these results indicate that these realizations are 
insensitive to temperature and power supply variations. A 
wide field of applications exists for the programmable filter 



beside 


the one discussed in 


Chapter VII. 


(1) 


Word recognition and 


speech synthesis; 


(2) 


Music applications; 




(3) 


Signal processing in 


communication ; 


(4) 


Adaptive balancing. 
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Further investigation is needed to develop a progranunable 
switched capacitor realization that can allow frequency 
scaling by changing clock frequency. Work can be also 
extended for developing a wide bandwidth progranunable filter 
using the composite operational amplifier technique proposed 
by [39]. Such implementation would lead to a very useful 
monolithic device at moderate cost. 

The research has yielded a paper that was presented at 
the 19th Annual Asilomar Conference on circuits, systems and 
computers, Monterey, California; November, 1985 (Appendix C). 
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APPENDIX A 






dimension ATlSf 1001 ,PT1S( lOO) ,AT2 3( 100) ,PT2S( I 00 I 
AT3S t lOOJ f PT35 f I 00 ) ,FAT l( 100 ) fFAT2 ( 1 00 J ,FAf3 i 100) 

O I ME NS I ON AT M lOO) y PT 1 I 100} >AT2O00 ) , PT2 t I 00 ) , A T3 C I 00 ) y 
^ PT3f lOO) ,FPTl( 100) yFPT2( 100) ,FPT3< 100) y 

^ FTFI (100) ,FTF2 (100) yFTF3 C I 00 ) y FP I ( 1 0 0 ) ,FP2 ( I OOl yFP3( 100) 
COMPLEX TIS yT2S yT33 yY 1 S y Y 2S y Y 3S y Y AS y Y 5S y Y6S y Y7S , Y 8S y D S , K 1 S , X2 S 
:^K3S,K ASyKSSy K6S 

COMPt.£X Tl yT2 yT3y Yl ,Y2| Y3y YAy Y5yY6y Y7y YByDyKl ,K2yK3 ,KA,K5,K6 
dimension FR( 100) y3TlS(l00),HTlS( 100) y3T 25(100) yHT2S( 100) 
<iy3T35( 100) ,HT3S( 100) yFBTK 100) yFBT2( 100) yF0T3( 100) 

D I PENSION BTl ( 100 ) , HT I ( 1 00} ,aX2( ICO ) , HT2 C I 00 ) y 9T3 i I 00 ) 

^yHT3f 100) fFHTli lOO) yFHT2( 100) ,FHT3C 100) 

COMPLEX CIS yC2S yC 3S , Gl S y G2S y G3 S y G AS , G5 S y G6S y G7S y G 0S , DOS y X 1 S yX2 
^X3Sy XAS y XSSy X6S 

COMPLEX Cl yC2 yC3y Gl ,G2yG3 yGAy G5yG6y G7yG8yODy XlyX2 yX3y XAyXSf X6, 
OMEGA = 0. 



DO 20 K = I y 100 

OMEGA = OMEGA ♦ 1 25 *^20 • 

S := CMPLX( 0.0 y OMEGA ) 

DATA 



PI = 1600. 

Cl = 200E-9 

RQl = 3200^ 

R2 = 1600. 

C2 = 50E-9 
RQ2 = 3200. 

R3 = 1600. 

C3 = 12.5E-9 

RQ3 = 3200. 

R4 = 1600. 

CA - 3.12E-9 
RQ4 = 3200. 

A3 = AA = A3S 
Al = A2 := AIS 
J = 4 

I = 4 
L = 4 
N = 4 



A4S = l£5 

A2S = lES 



^ ^ 4: :t: ^ ^ 









THE TYPE OF THE FOURTH FILTER 

DO CASE L 
CASE 
C HPF 

GIS = I ./R4 
G2S = GlS 
G3S = S^C4 
G4S = GIS 
G6S = GIS 
G5S = 0.0 
G7S = G3S 
G8S = l.yRQA 

CASE 

C LPF 

CIS = 1./R4 

G2S = S4:C4 

G3S = G2S> 1./RQ4 

G4S = GIS 

G63 = 0.0 

G5S = GIS 

GTS = 0.0 

G8S =GIS 

CASE 

C NOTCH 

GIS = 1./R4 

G2S = GIS 
G3S = S4C4 
G4S = GIS 
G5S = GIS 
G6S s 0.0 
GTS - G3S 
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G7S = 

Gas = 

I,- NONE OO 
8PF 

GIS = 

G2S = 

G3S = 

G4S = 

G5S = 

G6S = 

GTS = 

G8S =S^C4 
END CASE 



G3S 

1 ./ROA 



I •/RA 

G 2S 
S:>CA 
GIS 
0.0 
GIS 
I .yRGA 



SELECT THE TYPE OF THE THIRD FILTER 

OO CASE N 
CASE 

C HPF 



CASE 

LPF 



CASE 

NOTCH 



IF NONE 
8PF 



G1 


= 1./R3 


G2 


= Gl 


G3 


= S^C3 


GA 


= Gl 


G6 


= Gl 


G5 


= 0.0 


G7 


= G3 


G8 


= 1./RQ3 


Gl 


= 1./R3 


G2 


= S^C3 


G3 


= G2 ♦ X 


GA 


= Gl 


G6 


= 0.0 


G5 


= GX 


G7 


= 0.0 


Gd 


= GX 


1 

Gl 


= X./R3 


G2 


s Gl 


G3 


= S^C3 


GA 


= Gl 


G5 


= Gl 


G6 


= 0. 


G7 


= G3 


G6 


= X./RQ3 



X./RQ3 



OO 





GX 


= 1 ./R3 




G2 


= Gl 




G3 


= S^C3 




GA 


= G 1 




G5 


= 0.0 




G6 


= Gl 




G7 


= 1./RQ3 




G8 


= S^C3 


END 


CASE 




t:<t^«^^^j>«t^<r^SEL£CT THE 


DO CASE 


1 




CASE 

HPF 




YIS 


= 1./R2 




Y2S 


= YIS 




Y3S 


= SAC 2 




YAS 


= YIS 




Y6S 


= YIS 




Y5S 


= 0.0 




Y7S 


= Y3S 




Y8S 


= 1./RQ2 


CASE 

LPF 




YIS 


= 1./R2 
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Y2S = S^C2 

Y3S = Y2S^ 1./RQ2 

Y4S = YIS 

Y6S = 0*0 

Y5S = YIS 

Y7S 0.0 

Y8S =Y1S 

CASE 

C NOTCH 

YIS = 1*/’R2 

Y2S = YIS 
Y3S = S<^C2 
Y4S = YIS 
Y5S = YIS 
Y6S = 0*0 
Y7S = Y3S 
Y8S = 

IF NONE OO 
C 3PF 

YIS = 1 ./R2 

Y2S = YIS 
Y3S = S^2 
Y4S = YIS 
Y5S =0*0 
Y6S = YIS 
Y7S = 1 ./RQ2 

Y8S =S^C2 
END CASE 

SELECT THE TYPE OF THE FIRST FILTER 



OO CASE 


J 


CASE 

HPF 


Yl 


= l./Rl 


Y2 


= Y 1 


Y3 


= S<:Cl 


Y4 


= Yl 


Y6 


= Yl 


Y5 


= 0*0 


Y7 


= Y3 


Y8 


- 1«/RQ1 


CASE 

LPF 


Yl 


= 1«/R1 


Y2 


= S«C 1 


Y3 


= V2 ♦ l./RQl 


Y4 


= Y 1 


Y6 


= 0.0 


Y5 


= Yl 


Y7 


= 0.0 


Ya 


= Yl 


CASE 

NOTCH 


YX 


= 1 • /R 1 


Y2 


= Yl 


Y3 


= S^CX 


Y4 


= Yl 


Y5 


= Yl 


Y6 


= 0* 


Y7 


= Y3 


Y3 


= 1«/RQ1 


IF NONE 


DO 


SPF 


Yl 


= 1 ./Rl 


Y2 


= Y 1 


Y3 


= S*Cl 


Y4 


= Yl 


Y5 


= 0.0 


Y6 


= Yl 


Y7 


= l./RQ! 


Y0 


= S«C1 


END CASE 
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FR(K) = 0M£GA/6«2e 
S = CMPLXI 0.0, OMEGA) 

icio = ris*ti.» i./Aas) 

K2S = t Y2S •■ YSS ♦ Y6S) 
K3S = (Y4S ♦ Y7S ♦ Y8S) 
K45 = Y4S ♦ Y8S 
KSS = YIS ♦ Y3S 
K6S = YSS ♦ YSS 

K1 = 1./A2) 

K2 = |Y2 ♦• Y5 ♦ Y6) 

K3 = |Y4 ♦ Y7 ♦• Y8) 

K4 = V4 ♦ Y8 

K5 = Yl ♦ V3 

K6 = Y5 ♦• Y6 

XIS = G1S«(1.«’ I./A4S) 

X2S = <G2S *■ G5S ♦• G6S) 
X33 = tG4S ♦ G^S ♦ G8S) 
X43 = GAS ♦ G8S 
XSS = GtS ♦ G3S 
X6S = G5S ♦ G6S 

XI = GX«=( !.♦ l./AA) 

X2 = tG2 ♦ GS ♦ G6) 

X3 = 1 GA «• G7 ♦ G8) 

XA = G4 ♦ G8 

X5 = Gl ♦ G3 

X6 = GS •*• G6 



OS =<<1S«K2S«<3S)/A1S ♦ Y1S*Y4S*K6S ♦ Y3S**C 2S*K3S/ I A I S« A2 S) ♦ 

« Y2S^Y3S<«YAS/A2S«-( Y2S/'Y1S)PKIS#Y3S«( Y7S*Y8S) ♦ Y 3SSK6SSK3S/A 2S 

T1 S=( Y1S#YAS#YSS •• Y 3S* YSS« I Y7S- Y 83 ) *Y3S4tY7SP< Y2S^ YSS ) ♦ t Y7 5« 

# (Y1S»Y3S)«K2S)/AXSJ/DS 

T2S=iYlSPY5S*K4S - Y1 S^YYS# YSS ♦Y3S*Y5SS:K4SyA2S -Y3 S#Y6S«Y7 S/ A2 3 

* ♦ Y5SSY7S*K5S/Al S ♦ Y2SAY33PY7S ) /OS 

T3S= IK1S*Y7S*X2S/A1S ♦KISAYASAYSS ♦Y3S*Y7SCK2S/A 1SAA2S ♦Y3S4eYAS 

* *YSS/A2S ♦ I Y1S/Y2S)*K1S4Y3S*Y7S K6S* Y3S«J Y7S/A2S ) /OS 

O = ( K1 AK2^K3 1/ Al Yl»YA*K6 •*• Y 3*K2 *K3/ < A1 A2 ) ♦• Y2<:Y3#Y4/A2 ♦ 

« (Y2/Y1 >«Kl*Y3«l Y7«-Y8) •*• Y3«KS*K3/A2 

Tl = ( Y1 «YA#YS ♦ Y3AY54 ( Y7-Y8 ) ♦Y3#Y7A I Y2>Y6) Y7» | Yl Y3 ) «K2 ) / A 1 ) /O 

T2= 1 YX #Y5SsK4 -YX*Y7*Y6 ♦Y3«Y5#X4/A2 -Y 39Y64=Y7/A 2 ♦Y5«Y7«X5/A2 

# Y2*Y3^fY7 )/0 

T3=( KI «Y7«X2/A X ♦»< X#YA«YS ♦Y 3^:Y7«K2/AI* A2 ♦Y3*Y4«Y£/A2 ♦<Y2/YI)« 
A KXAY3AY7 ♦K6AY3AY7/A2 I/O 



OOS =( XXSAX2SAX3S)/A3S ♦ GXSAGASAXSS ♦G3SAX2SAX3S/I A3SA A4S) ♦ 

A G2SAG3SAGA3/AAS*(G2S/GXS) AXXSAG3SA(G73*G8S) ♦ G 3SAX6SAX3S/A4S 

Cl S= 4 GXSAG4SAGSS ♦ G3SAGSSA | G 7S-G 8S ) «-G3 SAG7SA J G2S^ G6S 7 ♦ I G7SA 

A (GXS^G3S)AX2S) /A3S) ZOOS 

C2S= ( GISAG5SAX4S -G1SAG7SAG6S •♦•G3SAGSSAX4 S/ AAS -G3SAG6S AG7 S/A4 S 
A ♦•G5SAG7SAX5S/A3S + G2SAG3SAG7S )/OOS 

C3S=tXXSAG7SAX2S/A3S ♦XXSAGASAG5S ♦G3SAG7SAX2S/A3SAA4S ♦G3SAG4S 
A AG5S/AAS ♦ 1 GIS/G2S7AXXSAG3SAG7S X 6SAG3SAG73/A4S )/OOS 

00 = 4 XIAX2AX3I/A3 GXAG4AX6 «-G3A X2AX3/ ( A3AA A ) ♦ G2AG3AG4/AA 
A XG2/G1 )AXl AG3A IG7^G8) ♦ G3AX6AX3/AA 

Cl>( GX AGAAGS + G JAG5A< G7-G8 )+G3AG7A( G2^G6) «-l G7A 4 G I + G3 ) AX2 ! /A3 ) /DO 
C2=( GX AG5AXA-G XAG7AG6 ♦G3AG5AXA/ AA -G3AGSAG7/A4 ♦G5AG7 AX5/ A4 
A G2AG3AG7 ) /OO 

C3 = t XX AG7AX2/A3 aXIAGAAGS ♦G 3AG7AX2/ A3AAA ♦G3AGAAG5/A4 ♦1G2/GDA 
A X1AG3AG7 ♦XSAG3AG7/AA i/OO 



ATISIK) = 20.AAX.0GX0ICABS( TIS) ] 
AT2SIKJ = 20 .AALaGIO<CABSI T2S) 1 
AT3S(K) = 20 .AAt_OGIO|CABS( T3S) ) 
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APTl S 
ART2S 
ART3S 

AI Tl S 
AI T2S 
AI T3S 



REAL (T15) 
REAL ( T2S ) 
REAL t T3S ) 

AI MAu( Tl SJ 
AI MAGI T2S) 
AI MAGI T3S) 



PTISK) = 
PT2SIK) - 
PT3SI>C) = 



ArAN2< AI Tl S » 
ATAN2IAIT2S t 
ATAN2IAIT3S • 



ARTIS)* 57.325 
ART2S)* 57.325 
ART3S)* 57.325 



ATlIKJ = 20.*ALaG10lCABSITl ) } 
AT2IKI = 20 .*ALOGl 0 lCAaS(T 2 )) 
AT3IK) = 20.*AL0G10ICA8SIT3) ) 



ARTI 

ART2 

ART3 

AITI 

AIT2 

AIT3 



REAL! Tl ) 

real! T2 ) 

REAL! T3 I 

AIMAGITI) 

AIMAGIT2) 

AIMAGIT3) 



PTinO = 
PT2IK) = 
PT3IK) = 

FATIIK) - 
FAT2IK) - 
FAT3IK) - 

FPTlIK) s 
FPT2IK) = 
FPT3IKJ = 



ATAN2IAITI 
ATAN2IA IT2 
ATAN2 I A ITJ 



, ARTI )* 57.325 

• ART2)* 57.325 

• ART3)* 57.325 



AT II K) 
AT2I K) 
AT3I K) 

PTH K) 
PT21 K) 
PT3I K» 



♦ ATIS(K) 

♦ AT2SIK) 

♦ AT3SIK) 



PTlStKI 
PT2SI K) 
PT3S IK) 



3TISHO = 20.*ALOGIO(CA8SIC13) ) 
BT2SIK) = 20 .PALOGIO ICA3SI C25) ) 
3T3SIK) = 20.*ALQG10(CABSIC3S) ) 

ARTIS = REALICIS) 

ART2S = REALIC2S) 

ART3S = REALIC3S) 



AI Tl S = 


AI MAGI Cl S) 








AIT2S = 


AIMAGt C2S) 








AIT3S - 


AI MAGI C3S) 








HTlSK ) 


= ATAN2I AI Tl S 


» 


ARTIS)* 


S7.325 


HT2SIK) 


= ATAN21 AI T2S 


1 


ART2S) * 


57.325 


MT3SIK) 


= ATAN21AIT3S 


i 


ART3S)* 


57.325 



3TUK) = 20. PALOGIOI CA8S1C I n 
8T2IK) = 20.CAL0G10ICA8SIC2) ) 
3T3IK) = 20.*AL0G10f CA8SIC3) ) 



ARTI 

ART2 

ART3 

AITI 
AI T2 
AIT3 



REAL! Cl ) 
REAL! C2 ) 
REALI C3 ) 

AIMAGICI) 

A1MAGIC2) 

AIHAGIC3) 



MTlIK) = 
HT2IK) = 
HT3IK) = 

FBTl IK) 
F3T2IK) 
FBT3 IK ) 



ATA M2 I AITI 
ATAIs.2 I A IT2 
ATAN2 I A IT3 



• ARTI )* 57. 325 

• ART2>* 57.325 

• ART3)* 57.325 



= 3THK) 
= 8T2IK) 
= BT3IK) 



♦ BTISIK) 

♦ BT2SIK) 

♦ BT3SIK) 
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c 



c 



c 



c 



c 

20 

60 



FHTICK) = HTUK) ♦ HTlSiK) 

FHT2IK) = HT2<K) ♦ HT2S(KJ 
FMT3tKj = HT3CK) ♦ HT3S(KI 

FTFXIX) = ATUK) ♦ STllK) ♦ BTlStK)* ATIS(K) 
FTF2tK) = FAT2XK1 ♦ FBT2tK) 

FTF3(<) = FAT34»0 ♦ F8T3(K) 



FPIIK) = FHTUK) ♦ FPTUKI 

FP2fK» = FHT24K) ♦ FPT2< K) 

FP34K) = FHT3IK) ♦ FPT31K) 



MRITEf 6«66) FR (K } ,FATU K ) ,F3T1 IK ) 

* ,FTFl nc) ,HTUtC) ,FHTl 

# FHT31K1 
CONTINUE 

F0RMAT<6| IX,F9.3J > 

STOP 

END 



GENTRY 
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APPENDIX B 



DIMENSION ATISI 10 0) fPTlSl 10 0) tAT2SC 100),PT2S( 100> 

, AT3S < I 00) , PT3S ( I 00 ) f FATU lOO ) ,FAT2 t 100) »FAT3 { I 00 ). 

DIMENSION ATl(l00)9PTl(l00)iAT2(100if PT2 (100)fAT3(l00)f 
^ PT3i 100) ,FPT1< lOO) ,FPT2< 100) ,FPT3l 100) » 

^ FTFl (100) f FTF2 I I 00 ) ,FTF3| 100 )♦ FPU 100) » FP2 I 1 0 0 ) t FP 3 ( I 00 ) 
COMPLEX TlS,T2SiT3SfYlS,Y2S fY3SfY45 f Y5S f Y 6S , Y7S , Y3S ,DSf KIS,K2S» 
^K3S,<4S|X5SfK6S 

COMPLEX r I tT2 ,T3f YX ,Y2f Y3, Y4» Y5 ♦ Y 6» Y7 » Y8 » D-r K I » K2 t K3 , K4 1 K5 f K6 
DIMENSION FR( lOOl , 3 T I S ( I 00 ) » H T I S { 100) ,8T2S{ 100) , HT2S C 100) 

<t^BT3S ( 100) ,HT3S| I 00 ) f F8 T 1 ( 1 00 ) , FB T2 ( I 0 0 ) »FBT3 ( 100 ) 
dimension OTl ( 1 00) t HTlt 100) tBT2( 100) tHT2( 100) tST34 100) 

^tHT3l lOO) tFMTU 100) ,FHT2C 100) tFMT3t 100) 
complex Cl S,C2S ,C3S f GlSfG2S>G3S,G4S ,G5S,<;65,G7S,G3S ,00S, XI S| X2S 
^?X3S fX4S f X5St X6S 

COMPLEX Cl^C2tC3fGl ,G2 t G3 t G4, G5,G6f G7 , G8 t DO , X I , X2 t X 3 , X4 ^ X5 . X6 , S 
OMEGA - 0. 

DO 20 K = 1,100 

OMEGA = OMEGA ♦ 125. ^20. 

5 2 CMPLXl 0.0 , OMEGA ) 

DATA 

Rl “ 1600. 

Cl - 2006-9 
RQI - 3200. 

R2 = 1600. 

C2 = 50E-9 
RQ2 = 3200. 

R3 = 1600. 

C3 = 12.56-9 
RQ3 = 3200. 

R4 = 1600. 

C4 = 3.125-9 
RQ4 = 3200. 

Wl = WIS = **S£T VALUE •• 

A3 - A4 = A3S s A4S = Vl/S 
Al = A2=A1S = A2S = WIS/S 
J - 4 
1=4 
L = 4 
N = 4 



C 

c 

c ^3j:^:C;3fc«:C:::M:^:C:^;>C^S£L£CT THE TYPE OF THE FOURTH FILTER 
DO CASE L 
CASE 

C HPF 

GIS = 1./R4 

G2S = GIS 
G3S = S<:C4 
G4S = GIS 
G6S = GIS 
G5S = 0.0 
GTS = G3S 
G8S = 1 ./RQ4 

CASE 

C LPF 

GIS = 1./R4 

G2S = S<^C4 

G3S = G2S^ 1./RQ4 

G4S = GIS 

G6S = 0.0 

G3S = GIS 

G7S = 0.0 

G8S =GIS 

CASE 

C NOTCH 

GIS = I./R4 

G2S = GIS 
G3S = S^C4 
G4S = GIS 
G5S = GIS 
G6S = 0.0 
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uu 



GTS - 
G8S = 

I.- NONE OO 
BPF 

GIS - 
G2S s 
G3S - 
G4S - 
G5S 5 
G6S = 

GTS = 

Gas =S^C4 
END CASE 



G3S 

1 



I ./RA 
GIS 
S^C4 
G IS 
0*0 
GIS 
I ./RQ4 



SELECT THE TYPE OF THE THIRD FILTER 

OO CASE N 
CASE 

: HPF 



CASE 

LPF 



CASE 

NOTCH 



BPF 



Gl 


= 1./R3 


G2 


Gl 


G3 


= S4C3 


G4 


= Gl 


G6 


= Gl 


G5 


= 0*0 


GT 


= G3 


G8 


= 1./RQ3 


Gt 


- 1./R3 


G2 


- S4C3 


G3 


= G2 ♦ 1 


G4 


= Gl 


G6 


= 0*0 


G5 


= Gl 


GT 


= 0.0 


G8 


- Gl 


1 

Gl 


= 1./R3 


G2 


= Gl 


G3 


= S^C3 


G4 


:= Gl 


G5 


= Gl 


G6 


= 0. 


GT 


- G3 


G8 


= 1./RQ3 


)NE 


OO 


Gl 


= 1 ./R3 


G2 


= Gl 


G3 


= S«C3 


G4 


=s Gl 


G5 


- 0.0 


G6 


= Gl 


GT 


:= 1./RQ3 


G8 


= S*C3 



. /RQ3 



END CASE 

c THE TYPE OF THE SECOND FILTER 



OO CASE 


1 






CASE 








HPF 


YIS 


— 


1 ./R2 




Y2S 




Y IS 




Y3S 


= 


S4C2 




Y4S 


— 


YIS 




Y6S 


s 


YIS 




Y5S 


s 


0.0 




YTS 


s: 


Y3S 




Y8S 


- 


1 ./RQ2 


CASE 

LPF 


YIS 




1 ./R2 
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CASE 

NOTCH 



Y2S 




S^C2 


Y3S 


= 


Y2S^ 


Y4S 




YIS 


V6S 


= 


0.0 


Y5S 




YIS 


Y7S 




0.0 


Y3S 


- 


YIS 


1 

YIS 


— 


1 ./R2 


Y2S 




Y IS 


Y3S 




S^C2 


Y4S 


= 


YIS 


Y5S 


s 


YIS 


Y6S 


s 


0«0 


Y7S 




Y3S 


Y8S 




1 ./R3 



1./RQ2 



IF NONE DO 
8PF 

YIS = 1 ./R2 
Y2S = YIS 
Y3S = s«:2 
y^S = YIS 
Y5S =0.0 
Y6S = YIS 
Y7S = 1./RQ2 

Y8S =S^C2 
END CASE 



SELECT 

DO CASE J 
CASE 

C HPF 



THE TYPE OF THE FIRST FILTER 



CASE 

LPF 



CASE 

NOTCH 



Yl 


= l./Rl 


Y2 


= Yl 


Y3 


= S<:C1 


Y4 


= Yl 


Y6 


= Yl 


Y5 


= 0.0 


Y7 


= Y3 


Y8 


= l»/ROl 


Yl 


= 1*/R1 


Y2 


= S<«Cl 


Y3 


= Y2 ♦ 1 


Y4 


= Y 1 


Y6 


= 0.0 


Y5 


= Yl 


Y7 


= 0.0 


Y8 


= Y 1 


1 

Yl 


= l./Rl 


Y2 


= Yl 


Y3 


= S<^C1 


Y4 


= Yl 


Y5 


= Yl 


Y6 


= 0. 


Y7 


= Y3 


Y3 


= l./RQl 



► /RQl 



IF NONE 
BPF 



DO 



ENO 



Yl 




1 ./Rl 


Y2 


= 


Yl 


Y3 




S4IC1 


Y4 




Y 1 


Y5 




o«o 


Y6 




Yl 


Y7 




l./RQl 


Y8 

CASE 


— 


S^Cl 
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c 

FRIK) = 0M£GA/6-2e 
S = CMPLXt 0.0, OMEGA) 

<13 = ris*ll.«- 1./A2S) 

<23 = ( Y2S ♦ Y5S ♦ Y6S) 

<33 = JY4S ♦ Y7S ♦ Y8S) 

<43 = Y4S ♦ Y8S 
<53 = YIS ♦ Y3S 
<63 = Y5S •• Y6S 
C 

<1 = Yl^=41.^ 1./A2) 

<2 = IY2 ♦ Y5 ♦ Y6) 

K3 = I Y4 ♦ Y7 ♦ Y8) 

<4 = V4 ♦ Y8 
K5 = VI ♦ V3 
<6 = V5 ♦ Y6 
C 

XIS S GlS4(t.*^ l./’A4S) 

X2S = <G2S ♦■ G5S ♦ G6S ) 

X33 - (GAS ♦ G?S ♦ G8S) 

X43 = GAS ♦ G8S 
X5S 3 GIS ♦ G3S 
X6S 3 GSS ♦ G6S 
C 

XI 3 G1<=(1.» 1./A4) 

X2 3 (G2 ♦ G5 ♦ G6( 

X3 3 (G4 ♦ G7 ♦ G8) 

X4 3 G4 ♦ GS 
X5 3 Gl ♦ G3 
X6 3 G5 ♦ G6 
C 

c 

03 3(K1S*K2S*<3S)/A1S ♦ YIS*Y4S*<63 Y3S3K2S6K3S/ ( A I S6 A2 S ) ♦ 

« Y2S4Y3S*Y4S/A2Sa( Y2S/’V1 3) *<1S#Y3S*I V7S>Y8S I ♦ Y3SS:<6S!><3S/A 2 
T1S=( Y1S4Y4S3YSS ♦ Y 3S* Y5S# { Y7S- Y 83 ) ♦Y3S#Y7S6( Y2S> V6S ) ■•■ ( Y7 33 

3 I Y 1SAY3S )3<2S) ✓AISJ /OS 

T2S3( YIS3V5S3<4S -YIS3Y7S3Y6S ♦ Y3S3Y5S3K4S/A2S -Y3 S3Y6S3Y7 S/ A2 3 

4 AYSSStVYSAKSS/Al 3 ♦ Y2S4Y334Y7S )/DS 

T3S3 (K1S4Y7S4K2S/AIS ♦■< I S4Y4 S4Y5 3 ♦ V3S*Y7S«<2S/A 1S4A2S ♦Y3S4Y43 
4 4Y5S/A2S ♦ I Y1S/Y2S)4K1S4V3S4Y7S ♦ < 634 V3S4 Y73/ A2S ) /OS 

C 

O 3| K14K24K3I/A1 ♦ Y14Y44<6 ♦ Y 34<2 4K 3/( A1 4 A2 ) ♦ Y24Y34Y4/A2 ♦ 

4 ( Y2/Y1 >4<1 4Y34 ( Y7*Y8) ♦ Y34K64K3/A2 

T13(Y1#ya 4Y5 ♦ Y34Yb4 1 Y7-Y8 ) *734 Y74 ( Y2* Y6 ) ♦ ( Y74 | Y I ♦ Y3 ) 4<2 ) /A 1) / 

T23(Y14Y54<4 -V14Y74Y6 ♦Y34Y54<4/A2 -Y34Y64Y7/A2 ♦YS4Y74K5/A2 «• 

4 Y24V34Y7 ) /O 

T3=( <l 4Y74<2/A 1 ♦KlAYAAYS aY 34Y74K2/A14 A2 ♦V34Y44YS/A2 ♦(Y2/Y1) 
4 <14V34Y7 ♦<64V34Y7/A2 »/0 

C 

c 

DOS 3( X1S4X2S4X3S) /A3S ♦ G1S4G4S4X6S >G3S4X2S4X3S/( A3S4 A4S ) ♦ 

4 G2S4G3S4G43/A4S*(G2S/G13) 4XlS4G3S4(G7SAGaS) ♦ G3S4X6S4X3S/A4 

Cl S3 ( G1S4G4S4GSS ♦ G3S4GSS4 ( G7S-G 8S ) >03 S4G7S4 ( G2 S«- G6S 7 ♦ ( G7S4 

4 (G1S«G3S}4X2S)/A3S)/OOS 

C2S3(G1S4G5S4X4S -G1S4G7S4G6S ♦G3S4GSS4X4S/A4S -G3S4G6S4G7 S/A4 S 
4 ♦G5S4G7S4X5S/A3S ♦ G2S4G3S4G7S )/OOS 

C3S3( X 1S4G7S4X2S/A3S ♦X1S4G4S4GSS ♦G3S4G7S4X2S/A3S4A4S ♦G3S4G4S 
4 4G5S/A4S ♦ ( G1S/G2S74X1S4G3S4G7S ♦ X 6S4G3S4G7S/ A4S )/OOS 

C 

003 ( XI 4X24X3 )/A3 ♦ G14G44X6 ♦G34 X24X3/ ( A34A4 ) ♦ G24G34G4/A4 ♦ 

4 (G2/G1 )4X1 4G34(G7«’G8) ♦ G34X64X3/A4 

Cl3( Gl 4G4 4 G5*G 34 G54 ( G7-G8 ) >034074 ( G2 AG6 » ♦ I G74f G 1 4G3 7 4X 2 7 /A 3 ) /DO 
C2=l Gl 4G54X4-G 14G74G6 >034054X4/ A4 -G34G64G7/A4 >G54G74X5/A4 > 

4 G24G34G7 ) /OO 

C33( XI 4G74X2/A3 >X14G44G5 >G34G74X2/ A34A4 >G34G44GS/A4 >(G2/G1) 

4 X14G34G7 >X64G34G7/A4 )/D0 

C 

c 

ATIS(K) 3 20.4A1_OG10(CA8S( TIS) 7 
AT2S(KJ 3 20 .4ALOG10(CABS( T2S) ) 

AT3SIO 3 20.4AL0G10(CABS( T3SI 7 
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c 



c 



c 



c 



c 



c 



c 



c 



c 



c 



c 



c 



c 



c 



c 



c 



c 



c 



APTIS = 
ART2S = 
ART3S = 


REAL (Tl SI 
REAL ( T2S I 
REAL ( T3S I 




AI TIS = 
AIT2S = 
AIT35 - 


A( MAG( Tl S) 
AI MAG( T2S) 
AIHAGC T3SI 




PT 1 S i < ) 
PT2S4K) 
PT3S(K) 


= ATAN2( AX Tl S 
= ATAN2(AIT2S 
= ATAN21AXT3S 


t ARTIS)* 57.325 
y ART2S)* 57.325 
y ART3SI* 57.325 


ATH K) 
AT21K1 
AT3t KJ 


= 20.#AL0G10lCAaS(Tl ) I 
= 20.«ALOG10(CAaS(T2n 
= 20««AL0G10<CA8S(T3} ) 


ARTl = 
ART2 = 
ART3 = 


REAL! Tl ) 
REAL! T2 I 
real! T3 I 




AITl = 
AIT2 = 
AIT3 = 


AmAG(Tl) 

AIMAG(T2) 

A1MAG(T3) 




PTl < K) 
PT2(K) 
PT3T KJ 


= ATAN2IAITI f 
= ATAN2(AIT2 t 
= ATAN2(AITJ « 


ARTl I* 57.325 
ART2I* 57.325 
ART3)* 57.325 


FATl iK} 
FAT2 (K) 
FAT3 CK) 


= ATUKl ♦ AT1S!KI 
= AT21 K) ♦ AT2S ( K) 
= AT3!K) ♦ AT3S(K| 


FPTl <K) 
FPT21 K J 
FPT3(<) 


= PT1!K) ♦ PTlStKl 
= PT2IKI ♦ PT2S(K1 
= PT3(K» ♦ PT3S!lC) 


arisiK) 

BT2S<K ) 
aX3S (K) 


= 20.:»ALOG10(CAaS! CIS} I 
= 20 .4SAL0GI0 (CAaS( C2S) 1 
= 20.4tALQG10(CAaS(C3S) ) 


APTIS = 
ART2S = 
ART3S = 


REAL! Cl SJ 
REAL (C2S I 
REAL IC3S) 




Al Tl S = 
AIT2S = 
AIT3S = 


AI MAG! Cl S) 
AIMAGI C2S) 
AIMAG! C3S1 




HTISIK ) 
HT2SK J 
MT3S4K) 


= ATAN2(AIT1S 
= ATAN2(AIT2S 
= ATAN2IAXT3S 


, ARTIS)* 57.325 
y ART2SI* 57.325 
y ART3SI* 57.325 


3TUK) 
BT2C K) 
aT3(K) 


= 20. «AL0G10f CAaSlC 1 } } 
= 20.*AL0Gl0!CAaS(C2) I 
= 20.AAL0G10!CA8S(C3) J 


APTl = 
ART2 = 
ART3 = 


real! C l 1 

REAL! C2 I 
REAL! C3 I 




ATTl = 
Al T2 = 
AIT3 s 


AIHAGXCll 

AIMAG(C2I 

AIMAG(C3) 




HTK K| 
MT21 K) 
HT3C K) 


= ATAN2(AIT1 ( 
= ATAN2(AIT2 t 
= ATAN2!AXT3 f 


ARTl 57* 325 
ART2)^ 57.325 
ART3J^ 57.325 


FBTl iKi 
FaT24K> 
F0T3 (K) 


= 3T1IK1 ♦ BTISIK} 
= 8T21K) «• BT2SOO 

= BT3!K» ♦ BT3S(K) 
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c 

FHTUK) = HTl(K) ♦ HTIS(K) 

FHT21KJ = HT2<K) ♦ HT2S(K} 

FHT3tK) = HT3CK) ♦ HT3S(K) 

C 

FTFllKJ = ATX(K) ♦ aTl(K) ♦ 6TISlK)> ATlS(Kl 
FTF2KI = FAT2IK) ♦ F9T2(KJ 
FTF3(KJ = FAT30CJ ♦ FBT3<K) 

C 

FPt(K) = FHTHK) ♦ FPTUK) 

FP2tK» = FHT24K) ♦ FPT2(K> 

FP3<K) = FHT3IK) ♦ rPT3(K) 

C 

WRITE (6,66) FR(K ) ,F AT U K) , F3T1 IK ) 

» ,FTF1I K) ,HTU K ) , FHTl (K) 

C ^ FHT3(K) 

20 CONTINUE 

66 FORMAT( 61 1X,F9. 3) ) 

STOP 

END 

GENTRY 
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DIGITALLY CONTROLLED "PROGRAMMABLE** ACTIVE FILTERS 



Sherif Michael and Panaglotla Andresakis 



Department of Electrical and Computer Engineering 
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Monterey, CA 9i9*J3 



ABSTRACT 

In this contribution a general purpose 
digitally controlled analog filter is presented. 
The novel design Is a cascade of second order 
sections that are Individually programmed to 
achieve any filtering topologies. Two-binary 
words are used to control the pole frequency wp 
and selectivity Qp of each section Independently. 
Each second order section is a Cenerallzed- 
Immlttance converter (GIC) blquads which are 
known for their high stability and low active 
and passive sensitivity. CMOS switches are used 
to electronically relocate the minimum number of 
passive elements to achieve function 
programmability. Switches are also used to 
select the number of cascaded sections to 
realize higher order transfer functions. 

1 . INTRODUCTION 

Tne availability of an analog filter with 
digitally controlled "programmable” coefficients 
has been the goal of many researchers due to Its 
several attractions. One possibility of a com- 
pact, versatile analog filter under remote 
control opens up many novel and independent 
application areas. Also, when a programmable 
filter is combined with a permanent reference 
memory which is user-progr ammabl e , this would 
form an economical and versatile device for 
dedicated stand-alone applications. The need 
for sticn a device was motivated by advancement 
In film and semiconductor technologies as well 
as the continuous upgrading of systems speclflca- 
tlona to taka advantaga of tna available 
technologies to tne limits. 

Linear analog filtering finds many appllca^ 
tlons, such as speech processing (recognition or 
synthesis), geology, instrumentation, coinmunlca- 
tlons, process control, adaptive balancing, etc. 
There has been much emphasis on performing the 
filtering function digitally, largely because of 
tne ease of varying and optimizing the transfer 
function. However, for many reasons, such as 
cost, size, signal processing complexity, and 
bandwidth. It would be desirable to perform the 
filter function with linear components yet 
retain tne flexablllty of varying the filter 
parameters digitally. 

Recently, tne advantages of combining 
linear components (operational amplifiers (OAS), 



resistors and capacitors) and nonlinear elements 
(switches) nave been demonstrated using switched 
capacitor techniques [1-3]. In this contrlcu- 
tlOM, we are presenting the results of realizing 
a continuous active device using linear elements 
and switches controlled by digital signals to 
achieve a fully programmable filter [^]. 
Several programming features of the proposed 
filter are reported. Tne first feature is the 
ability of the network to realize the most 
common filtering functions (function programma- 
bility) namely: Low Pass (LP), Band Pass (BP), 

Hlgi» Pass (HP), All Pass (AP) and Notch(N) 
functions, using the minimal set of elements. 
Tne second feature Is the ability of the network 
to program (Independently) tne key parameters of 
the filtering function chosen (parameter program- 
mability) namely: the pole resonent frequency 

(up) and selectivity (Qp). Finally the ability 
to program the network to cascade several sec- 
tions to achieve higher order filter. All of 
the above programmabl 1 1 ty features are performed 
Independently to realize a universal filtering 
network . 

2. DESIGN ANALYSIS OF THE PROPOSED FILTERS 

The basic active network considered. here as 
the heart of the programmable filter Is tne 
second order Generalized Immlttance Converter 
(GIC) structure [5]# Fig. 1, whose superior 
performance was established In tne literature 
[61. The general transfer function realized 
T(s) Is given by: 

T(s) •• M(3)/D(3)-(ao*a^3*as2)/(bo*bi3«b2s2) (i) 
Tne GIC transfer functions of Fig. 1 assuming 
non-ideal OAs are given by 




Stability and sensitivity analysis ; 



An Important criterion of a realization Is Its 
sensitivity to element variations. The GIC 
sensitivity analysis has shown to be as good or 
better than all competitive second order net- 
works [6]. While the GIC stability can easily 
be demonstrated since In all the transfer 
functions (2), the coefficients of D(s) are seen 
to remain positive for any OA mismatch. This Is 



. t 
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due to the absence of negative terms In D(s). 
Therefore the zeros of D(s) will remain in the 
left-half s-plane and low frequency unstable 
modes cannot arise during activation. 

Function programmability t 

The objective of this research was to 
develop a device that Is capable of realizing 
the following transfer functions: LP where T(s) 

• K/D(s), BP where T(s) - KS/D(s), HP where T(s) 
- KS2 /d(s), AP where T(s )-ls2-s{a)p/Qp) ]/D(s) 

and N where Ksl-Cs^^w^^l/DCs) . By optimizing 
the design of the filter, It was found that all 
of the above functions can be realized by the 
second order GIC section using four resistors, 
two capacitors and two OAs as shown In Table 1 . 
These passive elements are connected to differ- 
ent nodes to achieve the various realizations. 
A set of CMOS bilateral switches controlled by 
digital binary word, are used to relocate the 
same elements In different ways to achieve the 
desired filtering functions according to Fig, 2. 
The truth table of the switches control logic is 
shown In Table 2, while Fig. 3 Illustrates the 
corr espond Ing minimized CMOS logic circuit used 
for passive elements relocation. 

Paramater programmablll ty : 

While four of the resistors are equal and of 
value R each, the fifth resistor Is the Qp 
determining resistor and of value Rg-RQp. The 
two capacitors are equal and of value c-l/CupR) 
each. Two equal banks of capacitors are used to 
control*u)p. Each bank contains n binary- 
weighted capacitors connected In series with n 
CMOS switches as shown In Fig. . Using n bit 
binary word to control the switches, 2 ^ 
different values of o can be obtained that cor- 
responds to 2^ different values of up. Using a 
similar technique, the value of Rq can be con- 
trolled by an m bit digital word that yields 2™ 
different values of Qp, as illustrated In Fig. 5. 
Thus, full Independent control of the pole pair 
Up and Qp are achieved by programming the 
digital words controlling the switches to obtain 
the corresponding c and Rq. The complete second 
order programmable filter is shown in Flg.^.f, 
wher# the function programmability as wall as 
the parameter programmability are demonstrated. 

Higher order programmability ; 

Active filters design procedure can be 
classified as direct or cascade. In direct 
synthesis procedures the transfer function is 
realized as a single section [7]. In cascade 
synthesis procedures a high order transfer 
function Is expressed as a product of first and 
second order transfer functions and each of 
these is realized Independently. The overall 
network Is obtained by cascading the Individual 
sections. The cascade method of synthesis 
offers two practical advantages (a) simple 
network tuning (b) a few number of universal 
sections can be designed which can realize a 
multitude of network specifications. 



The second order GIC network structure 
lends itself to the cascade synthesis procedure 
since it does not require additional isolating 
amplifiers. Fig. 6b shows a block diagram of a 
programmable higher order filter that utilizes 
the second procedure by cascading 2 or more 
sections of the filter network shown In Fig. 6.a. 
The result is a high order fully programmable 
general purpose filter, that can be tailored to 
match almost any proposed specification, 

3. COMPUTER SIMULATIONS AND EXPERIMENTAL 
VERIFICATIONS; 



Fig. 7 shows dlfferenct computer simulation 
outputs of the programmable filter. The plots 
simulate the filter responses assuming ideal OAs 
with Infinite Gain Bandwidth Products (GBWP), as 
well as practical filter responses assuming OA‘s 
finite GBWPs of 1 MHz as of that of the LMT^I OA. 
A single pole OA model was utilized to approxi- 
mate the filter transfer functions In the later 
case. The approximation was found adequate 
since the simulation results of the nonideal 
response were found to be of close proximity to 
the experimental results of Fig. 8. The experi- 
mental results were obtained using a three bit 
word for filter topology programmability to 
select the type of transfer functions, A two 
words, four bits each, were used for filter para- 
meters programmability where wp and Qp are 
controlled Independently as given in Table 3. 
Fig. 7 also illustrates a higher order program- 
mability where a fourth order characteristics 
are shown for a LPF and a Chebychev BPF, 

^4. CONCLUSION 



The novel design described here has 
resulted In a universal programmable filter that 
can be digitally controlled to realize almost 
any practical filter specifications. This Is 
done through the use of CMOS switches controlled 
by binary codes to program the order of the 
filter, the filter topology, the filter center 
frequency and selectivity. The design procedure 
required developing optimum switching arrange- 
ments for the minimum redundancy in components 
and the least dependence of the filtering 
function on swltohlng Imp «r f «o 1 1 on a auoh aa 
switches stray capacitances and non-zero and 
non-llnear swltch-on resistance. Further invest- 
igation Is being conducted to develop a program- 
mable switched capacitor realization that can 
allow frequency scaling by changing clock 
frequency. Work Is also In progress for 
developing an extended bandwidth programmable 
filter using the composite operational amplifier 
technique proposed earlier by the author. Such 
Implementation would lead to a very useful 
monolithic device at moderate cost. 
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Fig. 4 Capacitors Banks for Frequency Programmability 
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Fig. 7 Computer Siimnulation' s Results. 
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Fig. 8 Experimental Results 
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